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ABSTRACT:  The  purpose  of  this  paper  Is  n  analyze  In  detail  many 
of  the  assumptions  and  techniques  often  employed  to  obtain  stress- 
strain  relations  at  rates  of  strain  of  approximately  lo3  sec"1. 

The  effects  of  axial  ware  propagation  and  Interaction  detail  In 
short  cylindrical  specimens  on  the  resulting  stress-strain 
relation  have  been  evaluated  In  detail.  Ihe  effect  of  specimen 
geometry  on  the  "end  effects*  associated  with  friction  has  been 
clearly  demonstrated  both  theoretically  and  experimentally.  An 
overall  evaluation  of  the  quality  of  the  data  generated  by  the 
split  Hopklnson  pressure  bar  experiment  based  on  a  one -dimensional 
wave  analysis  and  verified  by  experimental  observations  has  been 
presented. 
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INTRODUCTION 

The  behavior  of  materials  under  dynamic  loading  Is  receiving 
an  increasing  amount  of  attention  as  more  highly  transient 
environments  are  encountered.  In  order  to  simulate  the  rates  of 
loading  produced  by.impact  or  blast  wave  impingement, strain  rates 
from  50  sec~l  to  lCT*  sec“*  must  be  achieved.  Various  types  of 
impact  testing  devices  have  been  used  [1]*  which  not  only  permit 
but  require  the  existence  of  inertial  forces  cr  stress  wqvss  in 
the  specimen  in  order  to  determine  its  stress-strain-strain  rate 
relationship.  One  of  the  more  popular  devices  of  this  type  is  the 
split  Hopkinson  pressure  bar.  This  technique  consists  of  placing 
a  short  specimen  between  a  pair  of  bars  which  remain  elastic 
while  a  stress  wave  is  propagated  through  the  system.  When  one- 
dimensional  wave  analysis  is  applied  to  the  measured  strain 
records  from  each  bar,  both  the  force  and  velocity  at  the  specimen 
end  of  each  bar  can  be  calculated.  The  stress  and  strain  of  the 
specimen  at  any  lnstan~  of  time  is  obtained  by  averaging  the 
forces  and  velocities  obtained  from  the  elastic  bars.  The 
purpose  of  this  report  is  to  analyze  in  detail  many  of  the 
assumptions  which  must  be  made  in  order  to  draw  meaningful 
conclusions  from  this  technique. 

This  study  has  been  prompted  by  a  wide  variety  of  conflicting 
statements  which  have  recently  appeared  in  the  literature  which 
challenge  the  validity  of  the  split  Hopkinson  pressure  bar 
technique.  In  particular,  the  results  of  a  study  of  the  dynamic 
properties  of  a  high  purity  aluminum  utilizing  this  technique 
were  reported  by  Hauser,  Simmons  and  Dora  [ 2] .  They  concluded 
that  the  aluminum  tested  not  only  exhibited  a  dynamic  stress- 
strain  relation  but  that  the  relation  was  a  function  of  the  rate 
of  strain.  However,  various  experimental  analyses  [3,4]  of  wave 
propagation  in  this  material  indicate  that  the  rate  independent 
theory  of  wave  propagation  developed  by  von  Hainan  [53,  Taylor  [63 
and  Rakhmatulin  [7J  is  adequate  to  predict  its  response.  Bell  [4] 
has  recently  reported  an  experimental  study  of  the  split  Hopkinson 
pressure  bar  utilizing  his  diffraction  grating  technique  on  the 
specimen  to  obtain  an  independent  measure  of  strain.  Since  a 
large  difference  in  strain  was  observed  between  the  direct 
measurement  and  that  inferred  from  pressure  bar  measurements. 

Bell  concluded  that  "...the  source  of  the  difficulty  in  the 
extended  quasi-static  impact  tests  lies  in  the  assumption  of 
uniform  strain  in  the  short  elastically  bounded  specimen  and  in 
the  neglect  of  wave  propagation  and  interaction  detail." 

The  effects  of  the  wave  propagation  and  interaction  detail 
have  been  evaluated  numerically  in  this  report  by  applying  the 
method  of  characteristics  to  the  elastically  bounded  specimens. 


•numbers  in  brackets  designate  references  at  end  of  report 
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In  addition,  the  effects  of  friction  and  its  resulting  biaxial 
stress  condition  have  been  evaluated.  Various  experiments  were 
performed  at  a  variety  of  rates  of  strain  in  order  to  confirm 
certain  assumptions.  In  all  of  the  following  analyses,  the 
stress-strain  relation  of  a  material  will  be  assumed  to  be  unique 
in  order  to  simplify  the  calculations.  The  response  of  the 
pressure  bars  is  then  calculated.  In  order  to  obtain  the  stress- 
strain  relation  of  the  specimen,  the  same  assumptions  may  be 
employed  as  when  the  pressure  bar  response  is  obtained  experi¬ 
mentally.  The  usual  assumptions  are  as  follows: 

1.  Uniform  Axial  Stress  Distribution.  If  the  specimen  is 
made  sufficiently  short,  it  is  assumed  that  the  stress, and 
therefore  the  strain,  is  "effectively  the  same  throughout  the 
specimen."  [8]  This  assumption  is  equivalent  to  neglecting  the 
effect  of  axial  Inertia  in  the  specimen. 

2.  Prictionless  Interface.  The  presence  of  friction  at  the 
faces  of  the  specimen  will  cause  a  combined  stress  situation  to 
exist  at  the  interface  even  in  a  "static"  test  which  results  in 
an  "apparent"  increase  in  stress  for  a  given  strain. 

3.  Uniform  Radial  Stress  Distribution.  This  assumption 
effectively  requires  the  radial  stress  to  be  zero  and  the  axial 
stress  to  be  constant  across  the  diameter  of  the  specimen. 

The  Boundary  Conditions.  Consistent  with  a  one- 
dimensional  analysis  Is  the  assumption  that,  at  the  interfaces 
between  the  specimen  and  pressure  bars,  the  axial  forces  are 
equal. 

The  first  two  assumptions  will  be  analyzed  in  detail  in  this 
report.  The  remaining  assumptions  will  be  the  subject  of  a  two- 
dimensional  analysis  in  the  near  future. 


WAVE  PROPAGATION  AND  INTERACTION 

The  one-dimensional  "rate  Insensitive"  theory  of  finite 
amplitude  wave  propagation  has  been  utilized  to  determine  the 
transient  response  of  the  elastically  bounded  specimen  in  the 
split  Hopkinson  pressure  bar  experiment.  The  method  of  character¬ 
istics  has  been  used  to  obtain  a  solution  in  the  same  manner  as 
Conn  [93;  however,  the  solution  presented  here  has  been  carried 
to  a  logical  conclusion  which  permits  certain  generalizations  to 
be  made  regarding  the  effects  of  axial  Inertia.  The  approach 
taken  here  consists  of  assuming  the  stress  in  the  specimen  to  be 
a  known  unique  function  of  strain.  For  each  incident  strain 
pulse  the  wave  Interaction  detail  may  be  computed  at  any  point 
in  the  specimen, as  well  as  the  reflected  and  transmitted  strain 
pulsea  In  the  pressure  bars.  The  assumptions  outlined  in  the 
preceding  section  are  invoked  In  order  to  obtain  an  "apparent" 
stress-strain  relation  which  Includes  the  effects  of  axial 
inertia.  A  comparison  of  the  "apparent"  and  assumed  stress-strain 
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relations  indicates  the  error  to  be  expected  due  to  the  existence 
of  finite  amplitude  wave  propagation. 

The  governing  equations  which  describe  the  motion  of  the 
specimen  and  pressure  bars,  as  well  as  the  details  of  the  compu¬ 
tational  technique,  are  contained  in  Appendix  A.  The  calculation 
has  been  reduced  to  a  systematic  set  of  machine  instructions, 
the  CASH  code  (Appendix  B),  which  permits  the  automatic  computation 
of  the  entire  characteristic  net.  The  stress-strain  relation  of 
the  specimen  is  assumed  to  be  linearly  elastic  below,  and 
parabolic  above,  the  yield  stress.  However,  a  discontinuity  in 
the-  slope  of  the  curve  is  permitted  at  the  yield  stress  so  that 
a  wide  variety  of  materials  are  described  by  this  representation. 


Experimental  Confirmation 

Due  to  the  many  assumptions  contained  in  this  type  of 
analysis,  it  is  desirable  to  perfoim  z  calculation  where  experi¬ 
mental  data  are  available.  Therefore,  one  of  the  experiments  on 
commercially  pure  aluminum  recently  reported  by  Bell  [4],  which 
is  extremely  well  documented,  has  been  used  to  demonstrate  the 
accuracy  of  this  computational  technique  as  well  as  the  non- 
dimensional  nature  of  the  ons-dimensional  solution.  It  is  assumed 
that  the  specimen  is  linearly  elastic  to  a  stress  of  1000  psi  and 
the  resulting  stress  and  strain  at  yield  is  one  point  on  a 
parabola.  A  second  point  on  the  parabola  of  Bell  [4]  at  a  stress 
of  25,000  psi  has  been  used  to  obtain  an  adjusted  parabola  which 
when  normalized  with  respect  to  the  strain  at  yield  is  given  by: 

E  =  0.669*2  -  2.975207S  +  3.305787S2  (1) 


This  equation  is  presented  in  figure  1  where  it  is  compared  with 
the  parabolic  law  given  by  Bell  [4],  The  "hard"  aluminum  pressure 
bars  were  assumed  to  be  2024-T6  having  a  modulus  of  10'  psi  and  a 
density  of  J  lb/in3.  The  same  modulus  and  density  were  assumed 
for  the  specimen  in  the  elastic  region.  The  pressure  bars  were 
1.0  inch  in  diameter  and  the  specimen  had  an  t0/dp  of  \  and  a 
diameter  of  0.96  inch  which  closely  approximates  test  no.  827 
reported  by  Bell  [4].  The  data  deck  described  in  Appendix  B  for 
this  case  is  given  in  Table  1. 

The  computed  strains  at  the  center  of  the  specimen,  as  well 
as  the  computed  average  strains,  are  compared  in  figure  2  with 
experimental  values  reported  Xj  Bell  [4].  Relatively  good  agree¬ 
ment  is  obtained  between  the  computed  strain  and  the  average 
specimen  strain  obtained  from  pressure  bar  measurements.  Since 
dry  specimen -pres sure  bar  interfaces  were  used  for  this  test  and 
the  computation  assumes  perfectly  lubricated  interfaces,  the 
deviation  is  well  within  the  error  to  be  expected.  The  radical 
departure  of  the  diffraction  grating  measurements  at  approximately 
4  percent  strain  is  consistent  with  an  experimental  difficulty 
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reported  by  Bell  [10]  subsequent  to  the  publication  of  the  data 
contained  in  figure  2.  The  difficulty  is  associated  with  a 

Sronounced  change  in  reflectivity  of  the  surface  at  strains  above 
.3  percent. 


Axial  Inertia  Effects 

The  CASH  code  has  been  used  in  conjunction  with  the  REDUCE 
code  to  evaluate  the  effects  of  wave  propagation  in  the  specimen 
on  the  "apparent"  stress -strain  curve.  The  REDUCE  code  is  merely 
a  machine  program  which  is  used  to  process  the  incident,  reflected 
and  transmitted  strain  records  to  obtain  the  "apparent"  stress- 
strain  curve.  Although  originally  developed  to  process 
experimental  data,  it  is  directly  applicable  to  the  strain  pulses 
computed  by  the  CASH  oode. 

In  order  to  assess  the  various  effects  of  axial  Inertia  on 
the  resulting  "apparent"  stress-strain  curve,  the  results  of  the 
problem  considered  in  the  preceding  section  are  presented  in 
detail,  The  nonuniform  strain  distribution  which  is  most  severe 
during  the  earlier  portion  of  the  test  is  presented  in  figure  3(A) 
as  a  function  of  both  time  and  position  along  the  specimen.  In 
addition,  the  difference  in  stress  across  the  specimen,  which  is 
averageif  in  order  to  obtain  the  stress-strain  curve.  Is  presented 
in  figure  3(B).  It  may  be  seen  that  the  effect  of  nonlinear  wave 
propagation  in  the  specimen  is  to  damp  this  stress  difference 
which  oscillates  about  the  origin.  The  "apparent"  stress-strain 
curve  Is  presented  in  figure  4  with  the  assumed  stress-strain 
behavior  superimposed  for  comparison.  It  should  be  noted  that 
the  "apparent"  stress-strain  curve  Is  a  reasonable  approximation 
of  the  assumed  stress-strain  behavior  for  strains  greater  than 
1.5  percent.  The  "apparent"  overstress  occurring  at  strains  below 

I. 5  percent  is  an  axial  Inertia  effect  and  should  not  be  confused 
with  the  initial  overstress  which  occurs  during  the  Impact  of  a 
projectile  against  a  target  rod.  This  latter  effect  is  caused  by 
the  initial  three-dimensional  behavior  while  the  waveform  is 
established  in  the  target.  The  distortion  of  the  "apparent" 
stress-stain  curve  is  sufficient  to  preclude  the  determination 
of  the  yield  stress  or  modulus  for  the  case  considered. 

It  has  been  found  that  the  duration  of  this  initial  overstress 
is  approximately  twice  the  rise  time  of  the  incident  strain  pulse. 
Therefore,  great  care  must  be  exercised  in  analyzing  split 
Hopklnson  bar  data  obtained  during  the  time  required  to  establish 
the  initial  uniform  strain  distribution  in  the  sDecimen.  Tills 
problem  Is  alleviated  somewhat  If  the  material  being  tested  has 
a  well-defined  yield  stress  such  as  work-hardened  commercially 
pure  aluminum.  The  case  of  work  hardening  to  a  yield  stress  .f 

II, 000  psi  as  tested  by  Hauser,  et  al  [2],  and  analyzed  by  Conn  [ 9l 
has  been  considered.  The  assumed  and  computed  stress-strain 
curves  are  presented  in  figure  5.  Since  the  slope  oi  the  curve 

in  the  plastic  portion  of  this  material  is  much  less  than  in  the 
previously  considered  material,  relatively  large  nonuniform  strain 
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distributions  may  exist  in  the  specimen  with  a  much  reduced  effect 
on  the  stress  distributions.  Therefore,  a  good  approximation  of 
the  yield  stress  may  be  obtained  if  there  is  a  sharp  disconti¬ 
nuity  in  the  slope  of  the  curve  at  that  point. 

The  "acoustic  impedance"  of  a  rod  is  defined  as: 


I  -  00=A0  (2) 

It  has  been  found  that  the  ratio  of  the  pressure  bar  impedance 
to  the  specimen  impedance,  R,  is  a  sensitive  measure  of  the 
constancy  of  the  specimen  strain  rate.  The  impedance  ratio  was 
varied  by  using  steel,  titanium  and  aluminum  pressure  bars  in 
conjunction  with  various  specimen  geometries.  Figure  6  illus¬ 
trates  the  effect  of  this  ratio  on  the  average  specimen  strain 
rate.  The  wave  speed  required  to  compute  the  impedance  is  the 
plastic  wave  speed  which  is  a  function  of  time,  therefore,  the 
material  described  by  equation  (l)  was  assumed  In  all  'ases  and 
the  wave  speed  at  a  stress  of  10,000  psi  was  assumed  constant  to 
obtain  the  values  of  R  given  in  figure  6.  Since  this  ratio 
enters  the  calculation  as  a  single  quantity  in  the  boundary 
conditions  (Appendix  A),  it  is  reasonable  to  expect  that  a 
reduction  in  specimen  diameter  has  th~  same  effect  as  increasing 
the  density  of  the  pressure  bars.  This  is  reflected  in  figure  6 
by  the  results  for  ratios  of  15. ^  and  14.8.  The  former  are  the 
results  for  aluminum  pressure  bars  and  a  specimen  to  pressure  bar 
diameter  ratio  of  0.72  while  the  latter  are  the  result.*-  for 
titanium  bars  and  a  diameter  ratio  of  0.96. 

It  ha 8  been  shown  that  the  effect  of  averaging  the  end  face 
strains  will  yield  an  "apparent"  stress-strain  curve  which  is  a 
good  approximation  of  the  assumed  behavior  for  times  greater  than 
twice  the  rise  time  of  the  Incident  strain  pulse.  However,  If 
instead  of  averaging  the  pressure  bar  measurements,  the  stress  is 
assumed  to  be  constant,  then  the  transmitted  strain  pulse  Is 
proportional  to  the  average  specimen  stress  and  the  reflected 
strain  pulse  is  proportional  to  the  average  specimen  strain  rate. 
The  effect  of  axial  inertia  on  this  type  of  data  reduction  Is 
illustrated  in  figure  7  for  an  aluminum  specimen  described  by 
equation  (l)  with  an  t0/d0  of  unity.  The  oscillation  of  the 
curve  about  the  assumed  behavior  Is  exaggerated  because  of  the 
increased  length  of  the  specimen.  Therefore,  this  type  of  data 
reduction  is  only  useful  when  the  difference  In  the  stress  across 
the  specimen  Is  small  with  respect  to  the  average  stress.  Other¬ 
wise,  a  significant  error  could  be  Introduced  into  the  "apparent" 
stress-strain  behavior  if  an  attempt  is  made  to  average  the 
resulting  urve. 


MArnm  C.  ^  c  C 
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FRICTION  ANALYSIS 


The  analysis  contained  in  the  preceding  section  is  only  valid 
in  the  absence  of  frictional  forces  at  the  interfaces  between  the 
specimen  and  pressure  bars.  However,  compression  testing  of 
small  cylinders,  as  shown  in  figure  8(A),  at  any  speed  is 
extremely  difficult  since  the  combination  of  friction  and  specimen 
length  will  cause  a  nonuniform  stress  distribution  to  exist  near 
the  ends.  These  "end  effects"  will  give  rise  to  two  conditions 
which  are  experimentally  observable;  namely,  (a)  an  increase  in 
the  force  required  to  obtain  a  given  deflection,  and  (b)  a 
distortion  of  the  specimen  commonly  referred  to  as  "barreling." 

The  first  of  the  experimental  observables  previously  mentioned 
will  be  manifested  in  an  "apparent"  increase  in  stress  to  obtain 
the  same  strain.  The  following  approach  is  based  primarily  on 
the  analysis  reported  by  Jackson  and  Waxman  [ill.  The  magnitude 
of  this  effect  can  be  estimated  by  assuming  the  presence  of  a 
shear  stress  which  is  proportional  to  the  axial  stress  on  each 
face  of  the  specimen.  The  constant  which  expresses  the  ratio  of 
the  shearing  stress  to  the  axial  stress  is  defined  as  the 
"coefficient  of  friction"  and  is  given  by  the  symbol  u.  The 
differential  equation  which  expresses  the  static  equilibrium  of 
forces  shown  in  figure  8(B)  is  given  by: 


d57r 

3r“ 


+ 


(3) 


If  the  material  is  assumed  to  be  incompressible  in  the  plastic 
regime,  and  if  the  Treeca  yield  criteria  and  Hencky-MIses  flow 
law  are  employed,  it  can  be  3hown  [12]  that: 

=  ooe2“<s-r>A  (*) 

The  resulting  nonuniform  stress  distribution  giver,  by  equation  (4) 
is  shown  in  figure  8(C).  In  order  to  determine  the  effect  of  such 
a  stress  distribution  on  the  "apparent"  stress-strain  curve 
obtained  from  a  compression  test  with  friction,  the  force  must 
be  determined  as  a  function  of  specimen  geometry.  The  axial  force 
for  any  given  length  is  given  by: 

a 

c  rdr  (5) 

0  2 

If  the  "apparent"  stress  is  defined  as: 

aa  S  pz>a*  (6) 


p  =  2n  J 
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It  can  be  shewn  that: 


af/° a  “  Cea  -  (i-KX)] 


(7) 


where 


a  =  3  (1-Ht  )’3/2 


(8) 


Hare  3  la  the  ratio  of  the  coefficient  of  friction  to  the  original 
^o/^o  of  the  ap*ci*en  and  e  is  the  '‘engineering"  strain. 

The  increase  in  "apparent"  stress  as  a  function  of  strain 
given  by  equation  (7)  is  presented  in  figures  9  and  10  for 
various  values  of  the  parameter  3  of  practical  interest. 

Figure  9  demonstrates  the  error  to  be  expected  if  cylinders  with 
an  original  length  to  diameter  ratio  on  the  order  of  unity  are 
compressed  without  lubricating  the  anvils.  However,  even  when 
great  care  is  taken  to  maintain  a  frictionless  interface  with 
lubricants  such  as  molybdenum  disulfide  (u  =  0.04),  compression 
of  short  cylinders  with  an  (t0/d0)  «  .25  could  result  in  errors 
in  excess  of  5  percent. 

Associated  with  the  end  effects  of  a  short  cylinder  is  the 
distortion  of  the  specimen  commonly  referred  to  as  "barreling." 

This  is  a  two-dimensional  effect  which  could  cause  erroneous 
results  if  stresses  and  strains  are  deduced  from  force  and  cross¬ 
head  displacement  measurements.  However,  a  two-dimensional 
analysis  by  Davis  and  Jackson  [133  was  compared  with  the  preceding 
analysis  and  revealed  no  significant  differences  at  strains  below 
15  percent  [ 12] . 

In  order  to  confirm  the  validity  of  the  one-dimensional 
analysis,  several  experiments  were  performed  using  commercially 
pure  aluminum.  A  single  rod  of  1100-F  aluminum  as  received  was 
cut  Into  one-foot  lengths,  heated  to  650* F  for  2\  hours,  and 
oven  cooled.  All  specimens  used  In  this  study  were  machined  from 
this  one  piece  of  stock.  Cylindrical  compression  specimens  were 
machined  with  various  ratios  of  length  to  diameter  as  listed  in 
Table  2.  One  specimen  was  re -annealed  at  650*F  for  2\  hours  so 
that  the  effects  of  residual  stresses  due  to  machining  could  be 
evaluated.  The  faces  of  each  specimen  were  polished  with  emery 
cloth,  cesium  oxide  and  .Jewelers  rouge,  in  that  order.  A  pair  of 
compression  plates  of  4340  steel,  ground  to  a  number  three  finish, 
were  attached  to  a  standard  Baldwin-Southwark,  60,000-pound  capacity, 
universal  testing  machine.  The  crosshead  velocity  of  this  machine 
was  maintained  at  the  constant  values  indicated  in  Table  2  in 
order  to  maintain  the  same  constant  strain  rate  for  the  various 
geometries.  A  lubricant  consisting  of  molybdenum  disulfide  in  a 
light  oil  solution  was  used  throughout  this  series  of  tests. 
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Specimen  no.  1  was  loaded  to  a  atress  of  11,000  pal,  unloaded, 
re-lubricated,  and  then  reloaded  to  a  higher  stress.  This 
procedure  was  repeated  several  times  to  obtain  the  data  presented 
in  figure  11.  The  data  were  obtained  by  simultaneously  recording 
the  outputs  of  the  internal  balance  of  the  testing  machine  and  an 
electrical  deflectometer  attached  to  the  crosshead  of  the  machine 
which  was  assumed  rigid.  Specimen  no.  2  was  compressed  to  the 
same  final  strain  as  specimen  no.  1  without  re-lubricating  during 
the  test.  The  results  of  this  test  were  identical  to  those  illus¬ 
trated  in  figure  11,  indicating  that  lubrication  was  maintained 
throughout  the  test.  Specimen  no.  3  was  re-annealed  and  tested 
with  identical  results  indicating  the  lack  of  any  significant 
residual  stresses  due  to  machining.  Therefore,  all  remaining 
tests  were  performed  without  re-annealing  or  re-lubricating  the 
specimens  during  the  experiment.  The  remaining  specimens  were 
tested  and  the  "apparent"  stress-strain  curves  for  length  to 
diameter  ratios  from  2.70  to  0.277  were  obtained. 

The  force -deflection  data  collected  in  this  series  of 
experiments  are  presented  in  figure  12  for  three  representative 
ratios  of  length  to  diameter.  The  discontinuities  in  the  stress- 
strain  beh&Yior  were  repeatable  and  indicate  that  the  material 
considered  in  this  report  is  mechanically  unstable.  This 
observation  is  consistent  with  that  reported  by  Kenlg  [14]  where 
the  stress-strain  data  were  in  the  form  of  a  staircase  funotion 
instead  of  a  smooth  curve. 


In  order  to  determine  the  stress-strain  relation  for  a 
frictionless  test,  and  to  assess  the  accuracy  of  the  one- 
dimensiOnal  theory,  it  is  noted  that  the  exponential  in  equation  (7) 
may  be  represented  as  a  series: 


a 

e 


1 


+  a  + 


2  3 

a  a-5 

ST  +  3T  + 


(9) 


Substitution  of  this  expression  into  equation  (7)  and  neglecting 
higher  order  terms  yields: 


0/0 q  -  1  +  a/3 


(10) 


Therefore,  for  a  given  value  of  strain  the  apparent  increase  in 
stress  is  inversely  proportional  to  the  original  ratio  of  length 
to  diameter  of  the  specimen.  This  may  be  expressed  as: 


°1  -  °o  “l  81  <Vdo>2 

°2  -  °o  *  *2  "  *7  ‘  <Wl 


(11) 


Therefore,  the  stress  to  be  expeoted  from  a  friotionless  test  Is 
given  by: 
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<J  =* 

O 


(12) 


(t,/d0)0l  -  <Vdo>ff2 

<Wi  -  iw; 

Inserting  the  value 8  of  a,  and  for  t  /d  ’  s  of  2.7  and  0.5*5, 
respectively,  from  flgurei12  at  iO  pereent  strain,  it  is  found 
thatt 


(2-7)(1(^))-'(^:^{(17'05pl  »si  <»> 


or 


oo  -  15,080  pel  (1*) 

The  "apparent"  increase  in  stress  for  an  t Vd  of  .5*5  la  13.05 
pereent  at  a  strain  of  10  percent.  The  vaiue°of  P  which  will 
cause  this  increase  »ay  be  determined  from  figure  9  and  is  found 
to  be  equal  to  0.305.  The  coefficient  of  friction  may  then  be 
computed  from  the  definition  of  P  and  is  found  to  be  0.166.  The 
values  of  P  for  each  test  have  been  computed  for  this  coefficient 
of  friction  and  are  tabulated  in  Table  2.  A  computed  stress- 
strain  curve  for  p  equal  to  zero  may  be  calculated  from  the  data 
obtained  with  the  maximum  ^o/d0  (specimen  no.  6)  corrected 
according  to  equation  (7).  The  ^apparent"  stress-strain  curves 
predicted  by  the  one -dimensional  analysis  are  illustrated  in 
figure  12.  Excellent  agreement  is  obtained  with  the  Intermediate 
values  of  P  and  good  agreement  with  the  maximum  value  of  P 
considered.  The  maximum  deviation  of  the  computed  stress  from 
that  observed  experimentally  is  3.65  percent. 

These  results  indicate  that  the  effects  of  friction  and 
specimen  geometry  cannot  be  considered  independently,  but  will 
occur  in  a  predictable  manner.  In  addition,  the  measurement  of 
crosshead  displacement  is  an  acceptable  method  of  obtaining 
speoimen  strain  when  the  effects  of  friction  are  small  enough  to 
preclude  barreling.  Although  these  results  have  been  derived  for 
a  statloally  applied  load,  it  will  be  assumed  that  the  effects  of 
friction  and  specimen  geometry  are  synonymous  dynamically,  as 
well.  It  has  also  been  demonstrated  experimentally  that  even 
when  great  care  is  taken  to  eliminate  friction,  the  use  of 
relatively  short  specimens  will  introduce  an  "apparent"  increase 
in  stress  very  similar  to  reported  strain  rate  effects.  Therefore 
the  effects  of  friction  should  not  be  neglected  unless  it  can  be 
shown  that  the  geometry  of  the  specimen  can  be  changed  without 
affecting  the  resulting  stress-strain  curve  when  derived  from 
load  and  deflection  measurements. 
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EXPERIMENTAL  EVIDENCE  OP  DYNAMIC  BEHAVIOR 

If  the  material  under  investigation  is  sensitive  to  the  rate 
of  loading,  a  radically  different  response  may  be  indicated  if 
these  data  are  compared  with  the  data  obtained  from  low-speed 
testing  such  as  that  presented  in  the  preceding  section.  However, 
since  the  rates  of  strain  differ  by  approximately  six  orders  of 
magnitude,  such  differences  are  entirely  possible.  In  order  to 
demonstrate  that  this  different  response  is  a  property  cf  the 
material  and  not  due  to  the  presence  of  friction,  axial  or  radial 
inertia,  dynamic  data  were  obtained  for  commercially  pure  aluminum 
by  two  completely  Independent  test  techniques.  Medium  strain  rate 
data  were  obtained  with  the  more  conventional  high-speed  testing 
equipment  while  high  strain  rate  data  were  obtained  with  the 
split  Hopkinson  pressure  bar  technique.  Although  the  strain  rates 
produced  by  each  technique  do  not  overlap,  they  do  provide  data 
for  a  wide  variety  of  strain  rates. 


Medium  Rate  Testing 

The  first  series  of  experiments  to  be  described  was  tested 
with  the  high-speed  testing  equipment  shown  in  figure  13.  This 
apparatus  consists  of  a  modified  Plastechon  Model  581  testing 
machine  equipped  with  a  Schavitz-Bytrex  Model  PL  2500  load  cell, 
and  an  Optron  Model  68OAX  electro-optical  extensometer.  The 
extensometer  consists  of  two  optical  tracking  units  which  generate 
an  electrical  signal  proportional  to  the  displacement  of  the 
target  being  tracked.  The  targets  used  In  this  series  of  tests 
were  attached  to  the  upper  and  lower  faces  of  the  compression 
cage  shown  in  figure  13.  The  differential  output  of  the  two 
trackers  Is  a  voltage  proportional  to  the  "engineering"  strain 
as  a  function  of  time.  The  frequency  response  of  the  extensometer 
Is  reported  by  the  manufacturer  to  exceed  5  kc.  The  frequency 
response  of  the  load  cell -extension  arm-compression  cage 
combination  has  been  determined  experimentally  to  be  200  cps; 
however,  calibration  of  both  measuring  devices  was  accomplished 
statically.  The  output  of  both  the  load  cell  and  extensometer 
were  recorded  simultaneously  and  Independently  with  a  Tektronix 
Model  502A  dual  beam  oscilloscope.  The  data  taken  during  these 
tests  were  obtained  with  crosshead  velocities  sufficiently  low 
to  preclude  the  existence  of  axial  or  radial  inertia  effects  in 
the  specimens  and/or  "ringing"  of  the  measur_.ig  or  recording 
equipment. 

The  high-speed  testing  machine  operates  as  shown  schematically 
In  figure  14.  a  high  pressure  nitrogen  source  provides  a  constant 
pressure  to  a  piston  which  initially  compresses  the  oil.  By 
activating  a  solenoid  valve  the  fluid  is  allowed  to  discharge 
through  an  orifice  which  is  adjusted  to  obtain  various  constant 
crosshead  velocities.  The  compression  cage  is  fitted  with  a  slack 
adapter  which  permits  the  crosshead  to  accelerate  to  a  constant 
velocity  prior  to  engaging  the  specimen. 
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Three  specimens  were  machined  from  the  same  stock  used  for 
the  "static”  testing  program.  The  dimensions  of  each  specimen  were 
identical  (see  the  following  table),  and  the  same  lubricant  as 
used  in  the  "static"  tests  was  used  to  insure  that  any  differences 
in  observed  stress  could  not  be  attributed  to  frictional  effects. 
Hie  crosshead  velocity  of  each  test  was  varied  in  order  to  achieve 
the  variation  in  strain  rate  indicated  in  the  following  table. 

The  data  obtained  from  a  typical  test  are  presented  in  figure  15 
as  well  as  the  stress-strain  curves  for  the  three  rates  of  strain 
considered.  The  stress-strain  curve  derived  in  the  preceding 
section  for  equal  to  zero  also  is  illustrated  in  figure  15. 

The  data  obtained  from  specimen  no.  15  were  not  at  constant 
strain  due  tc  the  inability  of  this  type  of  equipment  to  maintain 
a  constant  velocity  at  high  rates  of  loading. 

Medium  Rate  Testing  -  Specimen  Configurations 
and  Test  Conditions 


Specimen 

no. 

Length 

(in) 

Diameter 

(in) 

Vdo 

Speed 
( in/sec ) 

e (sec-1) 

13 

0.500 

0.252 

1.98 

1.15-xlO"2 

.023 

14 

0.501 

0.251 

1.99 

8.75x10-2 

.175 

15 

0.50c 

0.251 

1.98 

0.5  to  1.5 

1.  to  3. 

Description  of  Split  Hopkinson  Bar  Apparatus 

The  general  arrangement  of  the  apparatus  used  in  this  series 
of  tests  is  similar  to  that  used  by  Krafft,  et  al  [ 153 ,  Llndholm  [lb], 
and  more  recently  by  Tanaka,  et  al  [17].  This  apparatus  is 
illustrated  schematically  in  figure  16.  The  pressure  bars  were 
mounted  on  four  Teflor  V-blocks  which  were  originally  aligned  by 
bore  sighting  from  the  high  pressure  section  of  the  projectile 
launcher.  Both  pressure  bars  were  machined  from  the  same  rod  to 
insure  the  same  properties  and  then  ground  to  the  same  diameter 
as  the  projectile  (0.483  inch).  The  material  used  was  7A1-4MO 
titanium  which  has  a  yield  stress  certified  to  exceed  172,000  psi. 

Each  pressure  bar  was  Instrumented  with  a  pair  of  Micro- 
Measurements  foil  strain  gages  (type  EA-06-125AD-±20)  diametri¬ 
cally  mounted  with  an  Eastrun  910  adhesive  for  strain  gage 
applications.  At  each  strain  measuring  station  the  gages  were 
wired  to  opposite  arms  of  a  Wheatstone  bridge  which  causes  the 
change  in  resistance  of  each  gage  due  to  axial  strain  to  be 
additive  and  that  due  to  bending  to  cancel.  The  output  of  these 
bridges  is  recorded  by  means  of  two  Tektronix  Model  531  oscilloscopes 
and  cameras.  The  incident  pressure  bar  was  instrumented  with  a 
Microsystems  semiconductor  strain  gage  (type  PA3-16-120)  located 
one  inch  from  the  Impact  point.  The  signal  from  this  gage  was 
used  to  trigger  two  delay  units  whicn  were  used  to  trigger  the 
recording  equipment  at  a  suitably  delayed  time.  Both  pressure 
bars  were  3b  inches  in  length  and  the  strain  measuring  station  on 
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eacn  pressure  bar  was  located  18  inches  froa  the  specimen  so  that 
the  simplified  method  of  data  reduction  introduced  by  Lindholm  [ 163 
could  be  utilized.  Integration  of  the  reflected  strain  pulse  was 
accomplished  by  means  of  a  Tektronix  Type  "0"  operational  amplifier 
plug-in  unit.  The  output  of  this  amplifier  was  used  to  drive  the 
horizontal  sweep  of  a  Tektronix  Model  531  oscilloscope  while  the 
transmitted  strain  signal  provided  the  simultaneous  vertical 
deflection.  The  calibration  of  the  Wheatstone  bridges,  opera¬ 
tional  amplifier  and  elastic  wave  speeds  is  described  in  detail 
in  Appendix  C. 

A  rectangular  stress  wave  is  generated  by  the  impact  of  a 
titanium  projectile  launched  from  a  small  bore  air  gun.  The 
projectile  used  was  10  inches  in  length  which  provided  a  pulse 
duration  of  approximately  100  microseconds  as  illustrated  in 
figure  17.  It  may  be  seen  that  a  rise  time  cf  approximately 
4  microseconds  is  possible  with  the  impact  of  a  flat  nosed 
projectile.  However,  the  resulting  oscillations  which  are 
introduced  due  to  the  finite  diameter  of  the  bar  are  undesirable 
for  several  reasons.  The  oscillation  of  the  incident  3train 
pulse  will  cause  an  oscillation  In  the  particle  velocity  at  the 
Incident  face  of  the  specimen.  This  condition  would  cause  a 
small  amount  of  nonuniform  work  hardening  on  one  side  of  the 
specimen.  Since  the  effect  of  cyclic  work  hardening  on  the 
dynamic  behavior  of  a  material  is  not  entirely  clear  [l6],  it 
should  be  avoided  if  at  all  possible.  In  addition,  the  averaging 
of  this  oscillation  is  only  valid  to  obtain  the  force  that 
generated  the  variation  in  Btrain  and  not  the  effect  at  the 
specimen  end  of  the  incident  pressure  bar. 

In  order  to  reduce  this  oscillation  as  much  as  possible  the 
impact  face  of  the  projectile  was  rounded  off  to  a  one-inch 
radius.  The  incident  strain  pulse  generated  by  this  type  of 
projectile  is  shown  in  figure  17(B)  for  a  slightly  lower  impact 
velocity  than  that  used  to  obtain  figure  17(A).  It  should  be 
noted  that  t’  e  rise  time  is  increased  to  approximately  10  micro¬ 
seconds;  however,  the  oscillation  is  attenuated  as  desired. 

Although  the  amplitude  of  this  incident  strain  variation  could 
be  attenuated  even  more  by  increasing  the  diameter  of  the 
pressure  bars,  the  variation  of  the  reflected  strain  pulse  would 
be  amplified  accordingly.  In  addition,  the  use  of  strain  gages 
with  relatively  long  gage  lengths,  which  has  been  suggested, 
would  only  serve  to  hide  this  variation  in  strain  and  not 
eliminate  it.  Neither  of  the  previously  mentioned  approaches 
were  attempted  to  further  reduce  this  oscillation. 

flpon  arrival  cf  the  incident  pulse  at  the  specimen,  part  of 
the  wave  is  reflected  back  into  the  Incident  bar  while  part  is 
transmitted.  The  transmitted  wave  propagates  into  the  "throw-off" 
bar  illustrated  in  figure  26,  reflects  from  the  free  end  and 
separates.  Just  as  the  "measuring  piece”  in  Hopkinson!o  original 
experiment.  This  technique  prevents  the  specimen  froa  being 
reloaded  by  any  additional  pulses,  thereby  permitting  its  recovery 
and  measurement. 
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The  specimens  used  in  this  series  of  experiments  were  again 
machined  from  the  same  stock  of  commercially  pure  aluminum 
previously  described.  Six  different  diameters  and  three  different 
lengths.  Table  3,  were  tested.  The  strain  records  obtained  with 
specimen  no.  19  are  considered  typical  and  are  shown  in  figure  16. 
The  same  lubricant  as  used  for  the  low- speed  tests  was  used 
throughout  this  program.  The  change  in  original  length  of  the 
specimen,  while  resulting  in  a  variety  of  strain  rates,  was 
Intended  to  demonstrate  the  lack  of  frictional  effects  tinder 
dynamic  conditions. 

As  described  by  Lindholrc  [ 16] ,  if  it  is  assumed  that  the 
stress  is  uniform  along  the  axis  of  the  specimen,  then  the  data 
may  be  reduced  electronically  as  shown  in  figure  16.  However, 
slnoe  this  involves  a  needless  additional  assumption,  the  data 
have  been  reduced  by  computing  the  forces  and  particle  velocities 
independently  at  each  face  of  the  specimen.  The  oscilloscope 
traces  were  "read"  on  a  Universal  Telereader  Type  17A  equipped 
with  a  Telecordex  which  automatically  punches  the  resulting  x-y 
coordinate  information  onto  IBM  cards.  The  three  data  decks 
(incident,  reflected,  and  transmitted)  are  used  as  input  to  the 
HSDUCE  code  which  numerically  determines  the  stress,  strain  rate 
and  strain  at  each  data  point.  This  program  also  computes  the 
difference  in  stress  across  the  specimen  in  addition  to  the 
final  stress-strain  curve.  The  output  of  this  program  for 
specimen  no.  19  is  presented  in  figure  18. 

Several  important  features  may  be  seen  by  comparing  the 
reduced  data  in  figure  18  with  the  raw  data  in  figure  1£>.  The 
relatively  small  oscillation  appearing  on  the  reflected  strain 
pulse  gives  rise  to  a  significant  oscillation  in  the  reduced  data. 
Since  the  frequency  of  this  oscillation  is  predictable  on  the 
basis  of  geometric  dispersion  in  the  pressure  bar  f 18] ,  this 
variation  in  strain  could  have  been  eliminated  prior  to  data 
reduotlen.  However,  the  alternate  method  of  retaining  this 
variation  in  raw  data  and  smoothing  the  reduced  data  was  utilized. 
This  problem  is  characteristic  of  the  type  of  material  being 
tasted  since  the  amplitude  of  such  oscillations  is  very  nearly 
proportional  to  the  amplitude  of  the  pulse  being  propagated  [ 19l . 
As  the  strength  of  the  specimen  is  increased,  the  amplitude  of 
the  reflected  wave  will  decrease  while  that  of  the  transmitted 
wave  increases. 

From  figure  18  the  average  strain  rate  is  seen  to  vary  from 
1450  to  900  «ee~*.  If  this  material  were  sensitive  to  the  rate 
of  strain,  this  variation  would  have  to  be  taken  into  account. 

0»e  technique  often  employed  is  to  cross  plot  stress  as  a 
function  of  strain  rate  for  each  value  of  strain.  This  requires 
several  tests  at  a  wide  variety  of  strain  rates  to  obtain  the 
stress-strain  curve  at  a  constant  strain  rate.  This  was  found  to 
be  unnecessary  for  the  commercially  pure  aluminum  tested. 

The  resulting  stress-strain  behavior  of  the  specimen  is 
illustrated  in  figure  18(C).  Based  on  ths  results  of  the 
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one-dimensional  ware  analysis  the  data  collected  during  the  first 
20  microseconds  (twice  the  rise  time  of  the  incident  pulse)  hss 
been  ignored.  The  oscillations  predictable  fro*  geometric 
dispersion  of  the  reflected  pulse  in  the  elastic  pressure  bar 
hare  been  smoothed  and  the  resulting  curve  superimposed  on 
figure  18(C)  for  comparison. 

In  order  to  demonstrate  that  the  increase  in  stx-ess  indicated 
by  these  experiments  is  independent  cf  the  specimen  length,  the 
stress-strain  beharlor  of  specimen  no.  21  is  compared  with  that 
of  specimen  no.  19  in  figure  19.  It  is  apparent  that  this 
material  exhibits  the  saiae  stress-strain  response  for  rates  of 
strain  fro*  1000  to  4000  iec~l.  However,  the  data  obtained  fro* 
this  series  of  experiments  are  relatively  scattered  and  cannot  be 
presented  as  a  single  curve  or  a  family  of  rate  sensitive  curves. 
The  variation  in  the  stress-strain  behavior  of  this  material  for 
this  series  of  experiments  is  illustrated  in  figure  19  and  ia 
Indicative  of  the  level  of  confidence  which  should  be  placed  in 
any  dynamic  results  for  this  material. 

Two  techniques  were  used  to  Independently  verify  the  accuracy 
of  these  experiments.  The  specimen  was  recovered  after  each  test, 
examined,  and  the  final  length  detenained.  Mo  evidence  of 
barreling  was  observed  and  the  final  strain  measurement  was  within 
one  percent  of  the  maximum  strain  indicated  by  the  resulting 
stress-strain  curve.  Although  this  indicates  that  the  specimen 
strain  may  be  computed  accurately  fro*  pressure  bar  measurements, 
it  does  not  provide  a  verification  of  the  stress  measurement. 
Therefore,  an  energy  balance  was  performed  in  an  attempt  to 
provide  this  verification.  Applying  the  theory  of  one-dimensional 
elastic  wave  propagation,  it  was  found  that  the  energy  absorbed 
by  the  specimen  is  given  by: 

®.  -  w-p  j;  <*i  -  -  ‘!>dt  <i5) 

where  r  is  the  duration  of  the  pulse.  The  energy  absorbed  by  the 
specimen  may  be  attributed  to  the  sum  of  the  strain  and  kinetic 
energies.  The  resulting  unbalance  in  energy  is  then  an  estimata 
of  the  accuracy  of  the  strain  measurements.  An  increase  in 
temperature  of  the  specimen  was  computed  assuming  an  adiabatic 
process  and  a  specific  heat  at  constant  volume  of  0.217  Btu/lb#P. 
Based  on  the  results  of  Farren  snd  Taylor  [20]  it  has  been  assumed 
that  93  percent  of  the  strain  energy  appears  as  a  temperature 
increase.  The  results  of  this  calculation  for  the  data  contained 
in  figure  16  are  listed  in  Table  4. 

Although  the  temperature  increase  of  6.28#F  indicates  that 
the  effect  of  thermodynamics  on  the  resulting  data  is  negligible 
for  this  particular  specimen,  it  does  not  provide  an  independent 
check  of  the  magnitude  of  the  stress-strain  z*elation  obtained 
f'om  this  test  since  the  same  balance  could  be  obtained  with  less 
strain  energy  and  only  a  slight  increase  in  internal  energy. 
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Table  4.  Energy  Analysis 

of  Specimen  No. 

19 

(in  lb) 

percent 

Incident  Strain  Energy 

122.85 

Incident  Kinetic  Energy: 

122.85 

Total  Incident  Energy: 

245.7 

Reflected  Strain  Energy: 

48.65 

19.8 

Reflected  Kinetic  Energy: 

48.65 

19.8 

Transmitted  Strain  Energy: 

15.5 

6.3 

Transmitted  Kinetic  Energy: 

15.5 

6.3 

Specimen  Strain  Energy: 

104.0 

42.3 

(AT  =  6.28#F) 

Specimen  Kinetic  Energy: 

0.06 

- 

Energy  Unbalance  Error) 

13.4 

5.5 

Total  Energy: 

245.7 

100.0 

18 
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DISCUSSION  OP  RESULTS  AND  CONCLUSIONS 

The  analyses  contained  In  this  report  have  been  performed  in 
an  effort  to  establish  the  limitations  of  the  split  Hopkinson 
pressure  bar  as  a  technique  to  obtain  dynamic  stress-strain 
relationships. 


Axial  Inertia  Effects 

In  the  analysis  of  results  obtained  from  the  split  Hopkinson 
pressure  bar  it  is  assumed  that  the  average  stress  and  strain 
rate  in  the  specimen  at  any  time  may  be  approximated  by  the 
average  of  the  stresses  and  particle  velocities  of  the  end  faces 
of  the  specimen.  The  effects  of  the  large  gradients  of  stress 
and  strain  caused  by  axial  wave  propagation  in  the  specimen  have 
been  evaluated.  Although  certain  conditions  have  been  found 
where  the  nonuniformities  of  stress  and  strain  are  severe,  the 
effect  of  the  averaging  process  is  to  yield  an  "apparent"  stress- 
strain  curve  which  is  a  reasonably  close  approximation  of  the 
actual  stress-strain  behavior.  Due  to  the  finite  rise  time 
required  to  develop  the  incident  strain  pulse,  a  significant 
overstress  will  be  apparent  during  the  earlier  portions  of  the 
test.  Fortunately,  this  overstress  is  not  an  accumulative  effect. 
However,  modulus  data  cannot  be  obtained  and  yield  stress  can 
only  be  obtained  approximately  by  this  technique  since  both 
phenomena  usually  occur  during  this  portion  of  the  test  when  the 
averaging  process  is  not  valid.  It  has  been  found  that  as  the 
ratio  of  the  "acoustic  impedance"  of  the  pressure  bar  to  that  of 
the  specimen  is  increased,  the  average  rate  of  strain  approaches 
a  constant  value.  The  effect  of  assuming  a  uniform  stress  and 
strain  rate  has  been  found  to  increase  the  amplitude  of  the 
oscillation  of  the  "apparent"  stress-strain  curve  about  the 
actual  curve.  However,  the  average  of  this  resulting  curve  agt.in 
appears  to  be  a  good  approximation  of  the  actual  behavior  of  the 
material . 


End  Effects 

The  effects  of  a  shearing  stress  on  the  faces  of  the  specimen 
have  been  examined  statically,  both  theoretically  and  experimentally. 
It  has  been  found  that  the  effects  of  friction  and  specimen 
geometry  cannot  be  considered  independently.  It  has  been  shown 
that  a  reduction  in  the  initial  to/d0  of  a  short  compression 
specimen  will  have  the  same  effect  as  an  increase  in  the 
coefficient  of  friction  for  the  same  specimen  geometry.  It  also 
has  been  demonstrated  that  "barreling"  will  be  initiated  at  lower 
values  of  strain  as  the  ratio  of  the  coefficient  of  friction  to 
the  specimen  ^(y/^0  is  increased. 

The  effects  of  friction  predicted  by  the  one-dimensional 
analysis  have  been  observed  experimentally.  Since  the  agreement 
between  the  predicted  load-defleotion  behavior  and  that  observed 
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experimentally  Is  considered  excellent,  a  method  Is  presented 
whereby  the  actual  or  frletionless  stress-strain  behavior  of  the 
material  may  be  computed. 


Experimental  Results 

Load-deflection  data  have  been  obtained  for  annealed  1200-F 
commercially  pure  aluminum  at  room  temperature  and  at  rates  of 
strain  from  1.8  x  10-3  sec-1  to  6.4  x  103  sec*1.  Low- speed  or 
"static"  testing  indicates  that  the  material  under  consideration 
is  rate  insensitive  at  the  lower  rates  of  strain.  In  addition, 
these  tests  indloate  that  commercially  pure  aluminum  is  a 
mechanically  unstable  material  as  defined  by  Kenig  [14], 

Reasonably  good  agreement  with  the  frictionless  low-speed 
test  results  was  obtained  With  a  medium  rate  testing  machine  at 
a  strain  rate  of  2.3  x  10“2  see”1.  However,  a  significant  change 
in  the  stress-strain  behavior  was  Indicated  as  the  rate  of  strain 
was  increased  to  .175  sec**1  for  the  same  specimen  geometry.  If 
the  results  obtained  from  the  preceding  friction  analysis  may  be 
assumed  valid  at  these  velocities,  then  the  Increase  in  stress  at 
a  given  strain  may  not  be  attributed  to  friction,  or  axial  or 
radial  inertia  effects.  It  was  observed  that  after  an  initial 
increase  in  stress,  the  data  approached  that  given  by  the  static 
response  of  this  material. 

Split  Hopkinson  pressure  bar  results  have  been  obtained  for 
this  material  at  a  variety  of  specimen  geometries.  It  has  been 
shawm  that  for  a  given  diameter,  the  length  or  the  specimen  may 
be  reduced  by  a  factor  of  four  without  significantly  affecting 
the  resulting  stress-strain  curve.  This  indicates  that  the 
effeots  of  friction  are  insignificant  at  these  rates  of  strain  if 
the  results  of  the  static  analysis  are  applicable. 

Due  to  the  scatter  in  data  obtained  for  this  material  at 
rates  of  strain  on  the  order  of  103  sec”1,  the  values  of  stress 
illustrated  in  figure  20  for  the  split  Hopkinson  pressure  bar  are 
only  considered  accurate  to  +5  percent.  This  relatively  large 
scatter  in  data  is  considered  to  be  indicative  of  a  definite 
limitation  of  this  test  technique  when  applied  to  mechanically 
cmstable  materials.  The  recently  reported  phenomena  of 
"catastrophic  straining  at  one  "point"  in  a  specimen,"  [21]  in 
annealed  1100  aluminum  could  easily  render  the  dynamic  testing  of 
short  specimens  meaningless.  This  phenomena  is  a  characteristic 
of  the  speolmen  and  is  not  indicative  of  the  accuracy  to  be 
expeoted  from  the  split  Hopkinson  pressure  bar  technique  with 
mechanically  stable  materials.  Probably  the  best  proof  of  the 
validity  of  this  technique  is  its  ability  to  determine  that  a 
stable  material  Is  insensitive  to  strain-rate;  such  as,  7075*-T6 
[22],  5he  dynamic  behavior  has  been  observed  to  be  the  sane  as 
the  static  relation  with  the  split  Hopkinson  pressure  bar  technique. 
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Finite  Amplitude  Wave  Propagation 

An  independent  method  of  obtaining  the  dynamic  stress-strain 
behavior  of  a  rate  insensitive  material  is  available  [23]. 
However*  before  considering  the  experimental  observations  it  is 
necessary  to  clarify  certain  misconceptions  which  have  developed 
over  the  past  few  years  regarding  the  "rate  independent"  theory 
of  plastic  wave  propagation. 

It  can  be  shown  [12]  that  the  same  differential  equations  as 
those  developed  by  von  Karman  are  applicable  to  a  rate  sensitive 
material  along  certain  paths,  i.e., 

dv  -  (16) 

poc 

along 


If  the  stress  is  expressed  as  a  function  of  strain  and  strain  rate 
then: 


v 


*  t|f).  +  (No 


de  c 


d  In  e 

<TT” 
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ir  th§  behavior  of  a  material  is  insensitive  to  the  rate  of  strain 
(da/oe  B3  0)  but  still  differs  from  the  "static"  behavior,  then 
the  measurement  of  finite  amplitude  wave  speeds  is  an  accurate 
measure  of  the  slope  of  dynamic  stress-strain  curve  at  each  level 
of  strain.  Therefore,  the  theory  employed  by  Bell  is  completely 
Justified  and  only  the  limitations  of  the  experimental  observations 
need  be  considered. 

Two  features  must  be  demonstrated  experimentally  in  order  to 
assert  that  a  material  is  insensitive  to  the  rate  of  strain.  First, 
it  must  be  shown  that  the  wave  speed  associated  with  each  level 
of  strain  is  a  function  of  the  strain  alone.  Second,  it  must  be 
shown  that  the  area  under  the  predetermined  wave  speed-strain 
diagram  will  uniquely  determine  the  particle  velocity.  The  first 
condition  may  be  demonstrated  by  relatively  few  experiments  since 
each  test  contains  strain  values  from  zero  to  the  maximum  strain 
determined  by  the  impact  velocity.  However,  the  second  condition 
requires  a  separate  test  at  each  velocity  for  which  the  rate 
insensitivity  is  to  be  demonstrated.  Unfortunately,  this  second 
condition  is  relatively  insensitive  to  changes  in  the  yield  stress 
of  the  material  because  of  the  violent  change  in  wave  speed  for 
very  small  values  of  strain.  Therefore,  in  order  to  demonstrate 
rate  insensitivity,  the  maximum  strains  developed  for  extremely 
low  impact  velocities  would  have  to  be  measured. 
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It  has  been  asserted  by  Bell  [24]  that  the  relationship 

c 


X 


sax 


c  dc 


(19) 


?ee/n  found  t0  apply  for  impact  velocities  "fro m  less  than 
100  in/sec  to  3000  in/sec. "  However,  since  the  amount  of 
scatter  in  the  data  that  had  to  be  averaged  to  obtain  these 
wave  speeds  does  not  appear  in  the  literature,  it  is  possible 
that  the  relationship: 


c 

max 

VD  +  J €  c  de  “  v0  (2°) 

D 

would  also  satisfy  the  same  experimental  observations.  If  this 
is  possible,  then  the  resulting  dynamic  stress-strain  relation 
could  be  found  from  the  expression: 


where  a  is  a  rate  sensitive  stress  occurring  at  a  strain,  c^, 
above  which  the  material  is  rate  insensitive.  Therefore,  the 
disagreement  In  the  dynamic  stress-strain  behavior  of  a  material 
obtained  by  the  measurement  of  wave  speeds  does  not  Invalidate 
the  results  obtained  from  split  Hopklnson  pressure  bar  experiments, 
but  does  demonstrate  the  rate  ineensitivity  of  the  material  for 
certain  strains.  These  results  are  in  agreement  with  recent 
experimental  observations  of  Bodner  and  Clifton  [25]  who  found 
the  plastic  defozmatlon  of  commercially  pure  aluminum  to  be 
Insensitive  while  the  yield  stress  was  sensitive  to  the  rate  of 
strain. 


CONCLUSIONS 

The  analyses  contained  in  this  report  permit  the  following 
conclusions  regarding  the  validity  of  the  split  Hopklnson  pressure 
bar  technique: 

1.  Experiments  utilizing  finite  amplitude  wave  propagation 

have  not  yet  demonstrated  the  Insensitivity  of  commercially  pure 
aluminum  at  strains  near  the  yield  point  and, therefore,  do  not 
contradict  various  split  Hopklnson  pressure  bar  results  appearing 
in  the  literature.  ^ 

2.  The  existence  of  nonlinear  wave  propagation  in  the 
specimen  significantly  affects  the  "apparent*  stress-strain 
behavior  of  the  specimen  for  a  time  equal  to  approximately  twice 
the  rise  time  of  the  incident  strain  pulse.  As  a  direct 
consequence,  it  is  unlikely  that  elastic  modulus  information 
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can  be  obtained  by  this  technique  for  most  metals.  In  addition, 
the  yield  stress  can  be  only  approximated,  unless  there  is  a 
distinct  change  in  the  slope  of  the  stress-strain  curve  at  the 
yield  stress. 

3.  "Apparent"  stress  and  strain  levels  obtained  at  later 
times  by  averaging  pressure  bar  measurements  are  not  signifi¬ 
cantly  affected  by  nonuniform  stress  and  strain  distributions. 

The  effects  of  friction  are  not  evidenced  in  the  "apparent" 
dynamic  behavior  observed. 

Commercially  pure  aluminum  has  exhibited  a  dynamic 
behavior  at  rates  of  strain  where  axial  and  radial  inertia  effects 
are  insignificant.  However,  split  Hopkinson  pressure  bar  results 
for  this  material  are  questionable  due  to  an  inordinate  amount  of 
scatter. 

It  is  concluded  that  this  technique  iur  the  determination  of 
dynamic  material  behavior  is  capable  of  generating  stress-strain 
relations  with  reasonable  accuracy  only  in  the  region  of  gross 
plastic  deformation  for  mechanically  stable  materials.  Maximum 
usefulness  in  terms  of  constancy  of  strain  rate  and  minimum 
oscillations  will  result  from  the  use  of  high  "acoustic  impedance" 
pressure  bars  with  a  minimum  diameter.  Although  the  effects  of 
friction  do  not  appear  to  be  as  severe  as  In  static  testing,  a 
lower  limit  on  the  specimen  t9/d0  must  exist  and  may  be  dictated 
more  by  radial  Inertia  and  adiabatic  heating  than  Interface 
friction.  This  is  a  problem  which  will  require  further  study 
but  does  not  affect  the  results  contained  In  this  report. 
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FIG .  1  STRESS-STRAIN  BEHAVIOR  OF  COMMERCIALLY  PURE  ALUMINUM 
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FIG .  3  EFFECTS  OF  AXIAL  INERTIA  ON  STRESS  AND  STRAIN  DISTRIBUTIONS 
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(A)  SPECIMEN  GEOMETRY 


(B)  DIFFERENTIAL  ELEMENT  OF  SPECIMEN 


(C)  STRESS  DISTRIBUTIONS  WITH  FRICTION 
FIG.  8  ONE-DIMENSIONAL  FRICTION  ANALYSIS 
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FIG.  9  ONE-DIMENSIONAL  FRICTION  EFFECTS  (0.1  </?<  0.5) 
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FIG.  10  ONE-DIMENSIONAL  FRICTION  EFFECTS  (0 </S< 0. 1) 
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FIG.  11  APPARENT  STATIC  STRESS-STRAIN  BE'  WIOR  OF  COMMERCIALLY  PURE  ALUMINUM 
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FIG.  12  EFFECTS  OF  FRICTION  ON  APPARENT  STRESS-STRAIN  BEHAVIOR 
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FIG.  14  SCHEMATIC  DIAGRAM  OF  MEDIUM  RATE  APPARATUS 


MEDIUM  RATES 


FIG.  16  SCHEMATIC  REPRESENTATION  OF  SPLIT  HOPKINSON  PRESSURE  BAR 
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(A)  IMPACT  OF  FLAT  NOSED  TITANIUM 

PROJECTILE  ON  TITANIUM  PRESSURE  BAR 


(B)  IMPACT  OF  SPHERICAL  NOSED  (1"  RADIUS) 
PROJECTILE  ON  TITANIUM  PRESSURE  BAR 


FIG.  17  EFFECT  OF  PROJECTILE  FACE  ON  INCIDENT  PULSE 
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FIG.  19  EFFECT  OF  SPECIMEN  GECMETRY  AND  SPECTRUM  OF  STRESS-STRAIN 
BEHAVIOR  AT  HIGH  STRAIN  RATES 
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FIG.  20  VARIOUS  OBSERVED  BEHAVIOR  OF  COMMERCIALLY  PURE  ALUMINUM 
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APPENDIX  A 

THE  ONE-DIMENSIONAL  ANALYSIS 

It  is  generally  conceded  that  the  one-dimensional  theory  of 
wave  propagation  will  accurately  predict  the  stresses,  strains 
and  particle  velocities  in  long,  thin  rods  where  the  area  change 
is  not  large.  However,  in  order  to  analyze  the  split  Hopkinson 
pressure  bar  experiment,  the  equations  must  be  applied  not  only 
to  the  pressure  bars  but  to  a  specimen  with  an  t0/d0  on  the  order 
of  unity.  Therefore,  the  response  of  the  specimen  must  be 
calculated  with  the  basic  assumption  that  a  uniaxial  stress 
condition  exists  at  all  times.  The  validity  of  this  assumption 
must  ultimately  be  established  experimentally. 


The  Pressure  Bar 

The  one-dimensional  analysis  when  applied  to  an  elastic  bar 
of  constant  original  cross-sectional  area  requires  that: 


dv  =  jJV  (A-l) 

O  o 

along  lines  in  the  a-t  plane  defined  by: 

a!  -  r  %  -  ±  Ca-2) 

Therefore,  along  those  lines  of  positive  slope,  integration  of 
equation  (A-l)  yields: 


v 


0  c 
o  o 


Similarly,  along  those  lines  of  negative  slope: 


(A-3) 


v  + 


c 


0  c 
o  o 


(A-4) 


The  constants  of  integration  K  and  K  may  be  determined  from  the 
Initial  conditions,  so  that  the  problem  is  reduced  to  the 
simultaneous  solution  of  two  linear  algebraic  equations. 

In  order  to  demonstrate  the  application  of  these  equations, 
consider  the  diagram  1"  figure  A-l.  It  is  assumed  that  all 
components  of  the  system  are  initially  at  rest  and  stress  free. 

An  arbitrary  pulse  has  been  generated  at  one  end  and  has  been 
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observed  at  the  strain-measuring  station  located  a  distance  Xj 
from  the  specimen.  The  constant  Kp  at  the  point  I  may  be 
determined  from  the  initial  conditions  and  is  equal  to  zero. 
Prom  equation  (A-3)  it  obviously  follows  that: 


vi "  °°€i 

o  o 


(A-5) 


The  value  of  K  may  be  determined  at  point  3  from  the  values  now- 
known  at  pointnI.  From  equation  ( A --4 )  it  is  seen  that: 


ctb  °t 

V  4-  5  =  V  4  — -  -  - 

B  po°o  1  °o°o 


(A-6) 


However,  from  equation  (A-5)  we  may  substitute  the  value  of 
into  equation  (A-6)  to  yield: 


v_  +  g-  =  2c  tT 
B  PC  o  I 

o  o 


(A-7) 


In  order  to  determine  the  value  of  Kn  at  point  R, it  is  necessary 
to  repeat  the  experiment  in  the  absence  of  a  specimen  with  the 
identical  incident  pulse  shape,  or  otherwise  generate  the 
incident  pulse  so  that  the  value  of  Kn  may  be  determined  in  the 
absence  of  reflections  from  the  specimen.  For  simplicity,  let  it 
be  assumed  that  Kn  may  be  determined  at  point  D  and  is  equal  to 
zero.  In  this  case. 


VR  “  '  °o€R 


(A-8) 


Applying  equation  (A-3)  between  points  B  and  R  yields: 


v  .  _!*_ .  v  _ 

B  Ooco  R  P0C0 


(A-9) 


Therefore, 


v  - —  =  -  2c  €_ 

B  PQc0  OR 


(A-10) 


Equations  (A-7)  and  (A-10)  may  be  solved  simultaneously  for  the 
stress  and  particle  velocity  in  terms  of  the  measured  quantities. 


-  'o<‘l  -  V 


(A-ll) 
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and 


°B  =  poCo2<eI  +  eR^  =  E^ei  +  e  R  ^  (A-12) 

The  Initial  conditions  of  the  problem  permit  equation  (A-3)  to 
be  applied  between  points  F  and  C  and  between  points  E  and  T. 
Therefore: 

(A-13) 

(A-14) 

Applying  equation  (A-4)  between  points  C  and  T  then  yields: 


v  -  _ — 

C  DC 


o  o 


VF  ~  D  c 


=  0 


0  o 


and 


aT 

VT  ~  p  c 
o  o 


"  VE  '  pc 


=  0 


o  o 


Vt 


(A-15) 


and 


\-/ 


p  c  2e 
00  T 


(A-l6) 


If  the  transmitted  strain,  e^,  is  measured  at  a  distance  Xm 
from  the  specimen,  then  the  stress  and  particle  velocity  or  the 
pressure  bar  at  the  specimen  at  time  tp  will  be  determined  from 

XT 

the  transmitted  strain  record  at  time  (t^  +  — ).  Similarly,  the 

C°  X-r 

pertinent  incident  and  reflected  strain  data  occur  at  (t0  - 

X  y  O 

and  (t_  +  — =-) ,  respectively.  Knowledge  of  the  stresses  and 
°  co 

particle  velocities  at  each  instant  of  time  then  permits  the 
stress  and  strain  in  the  specimen  to  be  computed  according  to  the 
definitions  contained  in  reference  8. 


The  Specimen 

Application  of  the  one-dimensional  theory  of  wave  propagation 
to  the  specimen  is  somewhat  more  complicated  since  the  paths 
along  which  the  differential  eouations  must  be  applied  are  no 
longer  straight  lines.  In  general,  the  differential  equations 
may  not  be  integrated  in  closed  form  and  the  many  wave  reflections 
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between  the  faces  of  the  pressure  bars  muse  be  taken  into  account. 
The  solution  to  this  problem  may  be  carried  out  graphically  and 
the  details  of  this  calculation  are  available  in  the  literature  [93. 
Therefore,  only  the  approach  to  this  problem  will  be  outlined 
here. 


The  stress  and  particle  velocity  at  an  arbitrary  point  in 
the  specimen  are  determined  by  the  simultaneous  solution  of  the 
two  differential  equations  which  apply  along  the  two  lines  that 
intersect  at  the  arbitrary  point  in  the  a-t  plane.  As  before: 


dv 


do 


poc 


however,  the  paths  are  now  defined  by  the  expression: 


(A-17) 


P 


2  de 

'  ST 


?o 

5T 


(A-18) 


For  a  rate  insensitive  material  or  for  a  constant  strain  rate 
test,  this  expression  may  be  simplified  to  yield  the  relation: 

o0c2  ,  (|f).  (A-19) 

Therefore,  along  the  positively  sloped  paths,  integration  of 
nuation  (A-17)  yields: 


v  -  v^  =  m 


(A-20) 


where 


dm  =  [p 


,(HV 


(A-21) 


Similarly,  along  the  negatively  sloped  paths,  integration  yields: 


v  -  v2  =*  -  <P  (A-22) 

At  the  boundary,  between  the  Incident  pressure  bar  and  the 
specimen,  the  magnitude  of  the  incident  stress  wave  in  addition 
to  the  waves  in  the  specimen  will  influence  the  motion  of  the 
interface.  At  any  point  on  this  boundary  the  characteristic 
paths  to  be  considered  are  a  positively  sloped  line  in  the  speci¬ 
men  and  a  negatively  sloped  line  in  the  pressure  bar.  In  the 
first  case  we  have  from  equation  (A-20): 
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VS  -  *5  =  V1 


CD. 


(A-23) 


while  in  the  second  case: 


B 


'B 


P  c 
00 


-  2Vi 


( (A-7) 


The  boundary  conditions  described  in  reference  8  may  be  expressed 
as  follows: 


vs  =  vB  (A-2l») 

and 


fs#s  *  °bab  <a-25> 

The  simultaneous  solution  of  equations  (A-7),  (A-23),  ( A-2^4 )  and 
(A-25)  yields  the  following  expression  from  which  the  specimen 
stress  may  be  determined: 

jlCf-  +  ®(®s)  =  2c0!z  -  Vj  *  (fi-26) 

Once  the  specimen  stress  is  determined,  the  velocity  may  be 
obtained  fi'om  equation  (A-23).  A  similar  treatment  may  be  applied 
to  the  specimen-transmitter  bar  interface  which  results  in  the 
following  expression: 


p-H-  *  a,(°S>  =  V2  -  ®2  'A-27> 

OOO 

which  may  be  solvea  for  the  specimen  stress.  This  stress  then 
may  be  used  to  determine  the  velocity  of  the  interface  from  the 
expression: 


vs  =  v2  +  -  cp(as)  (A-28) 

Providing  the  stress-strain  relation  for  the  specimen  and  the 
elastic  constants  of  the  pressure  bars  are  known  in  advance,  the 
stress,  strain  and  velocity  may  be  calculated  at  any  point  in  the 
a-t  plane,  including  the  boundaries.  Therefore,  knowledge  of  the 
incident  strain  pulse  permits  the  computation  of  both  the 


A-5 


NOLTR  67-156 


reflected  and  transmitted  strain  pulses.  Hiese  results  then  may 
be  used  to  compute  the  "apparent”  stress-strain  relation  which 
will  include  the  effects  of  nonlinear  wave  propagation  in  the 
specimen. 


The  CASH  Code 

Although  the  computation  of  the  stresses  and  velocities  in 
the  specimen  as  a  function  of  time  and  position  is  straightforward, 
the  graphical  method  used  by  Conn  [ 93  is  extremely  laborious  if 
any  degree  of  accuracy  is  desired.  Therefore,  this  problem  has 
been  programmed  in  FORTRAN  IV  so  that  the  task  of  evaluating  the 
effects  of  axial  wave  propagation  in  this  experiment  may  be 
lightened  somewhat. 

The  CASH  Code  (an  acronym  for  Characteristics  Applied  to  the 
Split  Hopkinson  Bar)  essentially  divides  the  a-t  plane  into  a 
number  of  lines  of  constant  a.  The  point  in  time  along  each  of 
these  lines  where  the  equations  are  to  be  evaluated  is  determined, 
as  shown  in  figure  A-2,  by  the  preceding  points  on  adjacent  lines. 
Hie  characteristic  lines  at  the  adjacent  points  are  extended 
assuming  no  change  in  slope.  In  general,  each  line  will  inter¬ 
sect  the  desired  line  at  a  different  time.  The  earliest  arrival 
establishej  the  location  of  the  point  and  the  other  line  is 
shifted  back  in  time  until  it  intersects  the  same  point.  Shifting 
of  the  adjacent  point  Is  accomplished  by  a  linear  interpolation  of 
the  function,  cp,  defined  by  equation  (A-21).  As  this  function 
changes,  the  slope  of  the  characteristic  line  also  changes  until 
the  resulting  line  intersects  the  desired  point.  The  particle 
velocity  then  Is  linearly  interpolated  to  the  same  point  and  all 
other  required  properties  then  are  evaluated. 

This  program  has  been  written  utilizing  several  subroutines 
so  that  other  problems  in  one-dimensional  wave  propagation  may  be 
solved  without  re -programming.  The  main  routine  is  used  to 
control  the  Input  of  data,  conversion  to  non-dimensional  quantities 
and  incrementing  through  the  characteristic  net.  Since  all  of 
the  equations  used  in  the  program  have  been  non-dimensionalized 
with  respect  to  the  specimen  yield  stress,  any  consistent  set  of 
units  may  be  used  for  input  data.  The  shape  of  the  Incident 
strain  pulse  may  be  described  by  as  many  as  200  pairs  (e,t)  of 
points  which  are  constant  or  increasing  in  magnitude  since  the 
present  subroutines  have  not  been  written  to  accommodate  elastic 
unloading.  Due  to  the  number  of  points  required  to  describe  the 
details  of  the  wave  phenomenon, only  ten  points  along  each  line 
of  constant  a  are  retained  at  any  time.  After  the  tenth  point  is 
computed,  all  ten  points  along  lines  for  which  output  is  desired 
then  are  stored  and  the  next  ten  points  may  be  computed. 

Subroutine  STORE  is  used  to  control  the  form  of  the  output 
of  data.  This  output  may  either  be  non-dim-nsional  or  have  the 
same  dimensions  as  the  input  data  at  the  option  of  the  user. 

After  ail  quantities  have  been  stored  at  50  points  along  as  many 
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as  eleven  lines  of  constant  a,  the  results  for  each  line  then 
are  printed  separately. 


Subroutine  BOUND  solves  the  characteristic  equations  at  the 
boundaries  of  the  specimen.  Equation  (A-26)  which  must  be  solved 
at  the  specimen-incident  bar  interface  may  be  written  as: 


S  +  P  -  2^j  “  +  pi  (A-26) 

where  the  non-dimensional  quantities  are  defined  in  the  list  of 
symbols.  Similarly,  equation  (A-27)  which  must  be  solved  at  tne 
specimen-transmitter  bar  interface  may  be  written  as: 


*tV 


S  +  P  -  u2  +  p2 


(A-27) 


Therefore,  both  equations  (A-26)  and  (A-27)  are  of  the  same  form; 


KjS  +  P  «  K2 


(a -29) 


where  K.  and  Kp  are  constants  which  must  be  determined  at  each 
interface.  Subroutine  BOUND  computes  the  value  of  each  of  these 
constants  at  each  time  and  then  uses  a  portion  of  subroutine  WAVE, 
which  contains  the  stress-strain  relation  of  the  specimen,  to 
solve  equation  (A-29). 

Subroutine  STRESS  solves  the  characteristic  equations  at  all 
interior  points  in  the  a-t  plane.  It  first  establishes  the  point 
in  time  at  which  the  equations  are  to  be  solved  and  controls  the 
shifting  of  the  characteristic  line  on  one  side.  The  equations 

to  be  solved  take  the  form: 


U  -  P  -  U2  -  ?1  (A-20) 

and 


U  +  P  -  U2  +  P2  (A-22) 

This  subroutine  uses  another  portion  of  subroutine  WAVE  in  order 
to  accomplish  the  shifting  and  establish  the  new  characteristic 
quantities  which  are  required  to  solve  equations  (A-20)  and  (A-2?). 

Subroutine  WAVE  is  the  only  portion  of  this  program  which 
requires  use  of  the  stress-strain  relation  of  the  material .  The 
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subroutine  reported  here  is  based  on  a  stress-strain  relation 
which  has  the  form: 


t  =  S 


for  S  <  1 


(A-30) 


and 


E  «  B2  +  B2S  +  B^S2  for  S  >  1  (A -31) 

This  equation  was  used  since  many  materials  may  be  represented 
with  a  stress-strain  relation  of  this  general  shape  [  26]  and  all 
derivatives  and  integrals  are  easily  obtained  in  closed  form. 

The  relationship  between  P  and  S  becomes: 

[3B^(P  -  1)  +  (BP  +  2B,)1'5]*  -  BP 
S  -  - 2 - * - ^ - -  (A-32) 


and 

(Bp  +  2B-S)1*5  -  (B0  +  2Bo ) 1  * ^ 

p  .  !  +  _1 - 3 - - 3 -  (*-33) 

The  slopes  of  the  characteristics  are  given  by: 

7S  =  +  (B2  +  2B3S)‘^  (A-3*0 


Thererore,  the  first  portion  ef  this  subroutine  is  used  to 
evaluate  the  stress,  strain  and  characteristic  slope  for  any  gjven 
value  of  P«  The  second  portion  of  this  subroutine  is  used  by 
subroutine  BOUND.  When  equation  (A-33)  is  substituted  into 
equation  (A-29)»  a  cubic  equation  in  3  results  which  may  be 
solved  for  its  one  real  root.  The  third  portion  of  subroutine 
WAVE  is  used  to  shift  the  required  characteristics  determined  by 
subroutine  STRESS.  The  methed  by  which  this  3hift  is  accomplished 
may  be  demonstrated  by  considering  figure  A-2.  Subroutine  STRESS 
is  used  to  determine  the  point  Tl  n  an<*  identifies  the  point 
T*L-1,N  which  needs  to  be  shifted  until: 


L-1,N 


tl»n  - 


(a— 35) 


The  slope  of  the  required  characteristic,  (5,  is  determined  by 
substitution  of  equation  (A-3*0  into  equation  ( A— 33 )  and  assuming 
that  P  varies  linearly  between  and  This  results 

in  a  single  equation  for  the  unknowA  time,  Tl_2,N*  which  may  then 
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be  evaluated.  This  time  then  Is  used  to  compute  the  values  of  P 
and  V  assuming  a  linear  variation  between  Tr_i  and  the  former 
Tl-i  w.  The  stress,  strain  and  characteristic  *  slope  then  may  be 
computed  from  equations  (A-32),  (A-33)  and  (A-34). 

The  solution  of  equation  (A-29)  and  the  shifting  of  charac¬ 
teristic  quantities  requires  the  solution  of  a  standard  cubic 
equation  with  one  real  root.  Since  this  must  be  accomplished  at 
almost  every  point  in  the  characteristic  net,  a  function  subroutine 
has  been  included  which  evaluates  the  one  root  without  the  use  of 
time-consuming  iterative  techniques. 

The  CASH  Code  provides  a  very  rapid  method  for  computing  the 
one-dimensional  response  of  the  eplit  Hopkinson  pressure  bar  sinoe 
there  are  no  iterative  schemes  employed  in  the  program.  It  is 
difficult  to  estimate  the  accuracy  of  this  program  because  of  the 
linear  interpolations  used  throughout;  however,  as  the  number  of 
points  in  the  characteristic  net  is  increased,  the  exact  solution 
should  be  approached.  The  sample  problem  described  in  this 
report  was  computed  with  MAX  equal  to  2C1  and  again  with  MAX 
equal  to  51  with  no  apparent  degradation  of  results.  Since  the 
computer  time  required  is  roughly  proportional  to  the  square  of 
MAX,  the  lower  value  is  recommended  for  reasons  of  economy.  Yne 
solution  to  the  sample  problem  was  carried  out  to  a  dimensionless 
time  of  60  with  MAX  equal  to  51  and  required  4.23  minutes  on  the 
IBM  7090  Version  13.  A  listing  of  the  FORTRAN  statements  whioh 
constitute  the  CASH  Code  is  provided  in  Appendix  B. 
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THE  CASH  CODE 


loFTC  CASH  LIST 

this  program  computes  Tr>t  transient  response  of  the  split 

HOPS I nSON  PRESSUre  EAR  EXPERIMENT  bY  THt  sfc THOU  OF  CHARACTER¬ 
ISTICS  ASSUMING  The  RADIAL  STR£SS  IS  EQUAL  To  7 FRO.  THF  MA I N 
ROUTINE,  CASH,  FORMS  ALL  VARIABLES  INTO  NON-L I  MENS IONAL  QUANTI¬ 
TIES.  ESTAml I ShEl  T m1  70N1NG.  AND  CONTP«AS  THE  INDEXING  THROUGH 
ThE  CHARACTERISTICS  NET.  The  FOLLOWING  DICTIONARY  OF  TERMS  MAY 
PROVE  HELPFUL  IN  MAXING  CHANGES  TO  THIS  PROGRAM. 

PROuR Art  Symbol  DEFINITION 

S  STRESS/SPECIhEN  YlcLL>  STfitSS 

U  PARTICLE  VELOC I TY/ImPaCT  VELOCITY  AT  YIELD 

E  ST.RAIN/SPrClMEN  YIELD  STRAIN 

C  WAV?  Vf LOC I TY/ELAST !C  WAVE  VELOCITY 

P  IMPACT  VELOC 11 Y/IMPACT  VELOCITY  at  YIELD 

T  IImE/TIMc  FOR  ELASTIC  wAVt  TO  TRAVERSE  SPECIMEN 

The.  FOLLOWING  INSTRUCTIONS  ARc  TO  ASSIST  Tnc  USER  I  i  PREPARING 
ThE  Data  DfcCX.  aNY  CONSISTENT  SET  OF  UNITS  MAY  E>E  USED... 

CI.E.  IN.  PSI,  IN/SEC.  ETC.) 

COL 

CARD  1  1-14  D 1  INCIDENT  BAR  DIAMETER 

15-2*  El  INCIDENT  BAR  modulus 

29-42  RhOI  INCIDENT  BAR  uEnSITY 

CARO  2  1  —  1 A  D2  TRANSMITTER  BAR  DIAMETER 

13-28  E2  TRANSMITTER  MODULU'' 

29-42  RH02  TRANSMITTER  far  DENSITY 

Caro  3  l-)4  DS  SPECIMEN  DIAMFTFR 

15-28  FS  SPECTMFN  MODULUS 

29-42  RHOS  SPECIMEN  DENSITY 

43-3o  SY  SPECIMEN  YIEL'u  STRESS  (MUST  HE  NuN-ZERu) 

CARO  it  1-14  G  DIMENSIONAL  CONVERSION  FACTOR  ( 386.0«6A  IN/ 

S£C**2 ,  OR  l.on  C'VSrC*»2) 

1S-28  XO  ORIGINAL  SPEC'mFN  LENGTH 

CARD  s  1-14  B 1  ARBITRARY  CONSTANTS  WHICH  DESCRIBE  THE 

1 S— 28  B2  STRESS-STRAIN  RELATION  OF  The  SPECIMEN 

29-42  83  IN  NUN-D!MENSiONAL  UMTS  (SEE  SUBROUTINE 

43-56  84  WAVE). 

F7-70  85 

CARD  6  1-72  ID  ANY  72  CHARACTERS  WILL  °F  PRINTED  FOR 

IDENTIFICATION  PURPOSES 

ANY  INTEGER  FROM  2  TO  11  WHICH  INDICATES 
ThE  huMoER  OF  POSITIvnS  ALOnG  THE  SPtCJMEN 

That  output  is  desired 

ANY  OUO  INTEGr.R  FROM  3  TO  201  WHICH 
DICTATES  THE  SIZE  OF  The  ChaP aC TF R I S T  I C  NET 
ENTER  1  IN  COL  15  IF  OUTPUT  IS  TO  BE 
DIMENSIONLESS.  leave  8LANX  JF  OUTPUT  IS 
DES1RFD  IN  SA«F  UNITS  AS  INPUT. 

CARD  8  1-5  X  ThE  ACTUAL  NUMBER  OF  PAIRS  OF  UATA  POINTS 

TO  BE  read  IN  (MAY  NUT  EXCEED  200) 

6-19  Cl  A  CONSTANT  WHICH  WILL  CONVERT  THE  INPUT 


CARD  7  4-5  I 

8-10  LMAX 
15  M 


B-l 


NOLTR  67-156 

APPENDIX  B 

DATA  TO  UNITS  OF  STRAIN 

20-33  C?  A  CONSTANT  WHICH  WILL  CONVFRT  THF  JNPUT 

DATA  TO  UNITS  OF  T I MF 

CAKU  913-16  £PS1(J)  4  PAIRS  OF  DATA  POINTS  PER  CARO  WHICH  ARE 

19-24  T!M(J)  PROPORTIONAL  TO  THE  STRAIN  AND  TIME  OF  THE 

INCIDENT  PULSE 

CARD  101-5  KKK  CONTROL  CARD  WHICH  DICTATES  THE  RETURN  FOR 

SUBSEOUFNT  CASES... 

-1  =  RETURN  TO  CARD  1 
0  *  RETURN  TO  CARD  3 
♦1  =  RETURN  TO  CARD  8 

COMMON  A,Al,A2.A3,A4,A5,ARLAl,ARbA2.B,Bl,B2,B3,B4,85,BE,C:201,12)» 
1C1  ,C2,C»AR1 ,CBAR2.CE,CO(  1 1 ,50  )  ,CR  ,CS  ,C'a  . D1  ,D2  ,DC  ,  DS  ,DT  ,DX  ,E(201 ,12 
21,61  »E2  »EF  ( 1  ?  1  ,F0<1 1 ,50)  ,FPS1  (20")  ,EPS21 12)  ,ER,FRO(  50)  ,FS,ETt  12)  , 
3ET0i  50 )  ,EW,G,I  .11 ,10(12)  .  J ,N , L  ,L'. ,LLL  ( 1 1  )  .LMAX.LMAX!  , L n.A X2,m»N»NN, 
4NMAX,nPAGE,NST0R£,P(201 ,12) ,P1 ,P2,P0( 1 1 , 50 1 ,PP ,PR .PW.RHCl ,RH02, 
5RH0S,S (20 1,121 .SLOPE, SMAX( 201 ) ,S0< 1 1 ,50 ) ,SR .S5.SW ,SY , T ( 201 , 12 ) . T1 , 
6  T2  »  T I  Mi  200 } ,T0( 11 ,50J ,TR,TS.U(201 .*2) ,U0( 11 ,50) .UR.US.XBARtl 1 ) ,X0 
200  READ  ( 5 , 1 i  D1,£1,RH01 
READ  (5,i)  D2.E2.RHC2 
205  REAO  (5.1)  DS.FS.RHOS.SY 
READ  (5.1)  G , XO 
READ  (5,1!  B1 ,B2,B3,P4,B5 
READ  (5.5)  ID 
READ  (5.2)  1,  LMAX.M 
210  READ  (5,3)  K.C1.C2 

READ  (5,4)  (EPS1 ( J) ,T!M(J) ,  J=1,K) 

WRITE  (6,6)  ID 

WRITE  (6.7)  D1,E1.RH01,D2.E2,RH02.DS.ES.RH0S,SY 

WRITE  (6,8)  B1 »m2.83»H4,H5 

WRITE  (6,9)  l.LMAX.K 

NSTORE  *  O 

NPAGE  «=  0 

DO  10  L *  1 »LMAX 

U( L , 1 ) =0, 

S( L  ,  1 ) *0, 

C ( L  » 1 ) *1 • 

E1L.1 1*0. 

T(L.l 1*0. 

P(L«1 1*0. 

SMAX ( L 1 *1 • 

10  CONTINUE 
EF ( 1 ) *0. 

ET< 11*0. 

EPS2 ( 1 ) *0. 

ARFA 1  * ( DS/D1 )«*2. 

AREA2* ( DS/D2 ) . 

C0ARli'SORT(Fl*RHOS/(ES*RHOl  )  ) 

CBAR2«3QRT(E2*RH0S/(ES*RH02) ) 

CS*SQRT ( ES*G/RHOS ) 

Al-SY/ES 
A2*CS*A1 
A3*X0/C  S 

A4*1.+RH0S*ARFA1/ (Rh02»C°AR1  ) 

A5-1 .+RHOS*AREA?/ (RH0?*CRAR2) 

C2*C?/A3 

LMAX1  *  LMAX  -  1 
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A  a  LMAX1 
LMAX2=LMAX1-1 
DX«1 ./A 
XBARI 11=0. 

1 1 =  1  —  1 

nn  jo.d 
8=FL0AT(J1/FL0AT(  1 1  1 
XBAR ( J. 1 1=8 
20  CONTINUE 

DO  350  11=1*1 

ULtlll  =  XBAR  (II)  »  FLOAT(LMAXl)  +  1. 

350  CONTINUE 
J=? 

SL0PE  =  C1*FPS1  (2)/(C2*dV(2l  1 

T(  1,?)=?.*DX 

EPS?! 2 l=SLOPF*DX*2. 

DO  60  L  =  2  .LMAX1 *2 
T ( L  *  1 ) =DX 
60  CONTINUE 
L  =  1 
N  =  ? 

S(1 ,21  =  (2.*0PARi*rp«;?(?)*Fe,}/{eY#()i+A?  fa  i  *Ri40<;/ (  tpaRi  *R«0]  I  1  1 
P(1.2l*S(1,21 
C ( 1 • 2 1  =1. 

E(l. 21  =  SI1.2) 

Ud  >2  )  =P d  *2  > 

IF  (Sd,?)-1, )  6A.6A.62 
62  SMAX(1 J=S(1.21 
6A  DO  70  l.  =  3*LMAX,2 
P ( L .2 1 *0. 

S  (  L ,2 1 =0. 

U(  L .2 1 =0» 

E ( L . 2 1 =0 . 

C ( L .2  1  =  1 . 

7  (  L  »  2  1  =*2  • *DX 
70  CONTINUE 

DO  BO  L=?,LMAX1  ,2 
P ( L . 2 1 =0, 

S( L  *2 1 =0. 

U ( L  *  2 1 *0. 

E ( L .2 1 =0. 

C  ( L  *  2  1  =  1  • 

TIL. 21=3. *DX 

80  CONTINUE 
P(2.21=Pd,21 
S  (  2  »  2 1 =S( 1.21 
U  (  2 , 2 1 =U ( 1 ,2  1 
F  (  2  *2 1 =E ( 1 ,2 1 

81  DO  110  N  =  3  *  1 2 
N  =  N 

L  *  1 

CALL  BOUND 

DO  90  L=3.LMAX2,2 

L*L 

NN=N-j 
CALL  STRESS 

90  CONTINUE 

91  DO  95  L=2.LMAX1,2 

IF  ( S ( L  *N— 1 1  -  SMAX(LI)  95,95.92 
9?  SB AX (LI  *  SIL.N-1 ) 
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9  5  CONTINUE 
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call  rouno 

DO  100  L'2  »l.MAX  1  i  2 

L  *L 

NN*N 

CALL  STRESS 

100  CONTINUE 

101  00  10*  L  = 1 ,LMAX,7 

IF  ( S(L .Kj-SMAXCL » )  10^,105,102 

102  SM4X<L)“S<L.N) 

105  CONTINUE 

IF  (TU.NI  .LT.  0.)  00  TO  990 
IF  ( S ( 1 «N )  ,LT.  0,1  GO  TO  980 
IF  (S(LMAX,N)  ,LT,  O.J  GO  TO  985 
IF  (J-K)  110, i lo.ooo 
110  CONTINUE 
CALL  STORE 
no  oc  no  l  =  i  ,lmax 
T(L»1»*T(L»11  J 
S(L  ♦n*S(L.ll  > 

UIL.l  I»UIL»n  > 

EtL»!)=ElL»ll I 
C(L*1!*C(L»I1) 

P.L.I )«P(LtlI> 

P(L.2)*P(t,12l 

T!L«2!*T!L,12I 

S(L»2)«S(L,12) 

U(L .2J=U«L,12I 
E(L»,)3E(L»I2l 
CIL.2)=C(L»12) 

1A0  CONTINUE 

FPs?m«cPS2n  1 1 
FPS2(2)*FPS2(12) 

FFI1  Infill  ) 

FF( 2 ) «FF ( 1 2  ) 

ET 1 1 )*ETI1 1 ) 

ET(21«ET(12I 
GO  TO  81 

980  CALL  STORE 
CALL  OUTPUT 
WRITE  t  6,9P1 I 
GO  TO  150 

985  CALL  STORE 
CALL  OJTPuT 
WRITE  (6,986) 

GO  TO  150 

900  CALL  STORE 
CALL  OUTPUT 
GO  TO  15n 

1  FORMAT  (5E14.51 

2  FORMAT  ( 3  >  5 1 

3  FORMAT  { I 5.2E3  4, 5 ) 

4  FORMAT  (12X.8F6.0) 

5  FORMAT  (12A6J 

'  -ORMAT  (1H1,12A6) 

7  FORMAT  ( 1H-,24X,RMDI AMETFR,12X,7HMC0ULUS,1 3X ,7HDFNSITY,1 1 X ,1 2HY I  EL 
10  STRESS/I 3uf|l NC IOFnt  BAR,3X,3r2n,8/I 6H0TR ansm j  TTPR  BAR,3E2o,8/ 
29H0SPFCIMrN,7X,4r?0.8! 

8  FORMAT  ( ’ H- , I  2 X , 2H3 1 , 1 8X , ?HB? , ; 8 X , 2HP 5 , 1 8X ,?HR4 , 1 8X  , 2HB5/5E20.8 ) 

9  FORMAT  (1H-.3HI  =,I5,20X,6HL«AX  *,I5,?0X,3HX  «,I5) 

981  FORMAT  ( 1H0,62H5PFCIMFN  HA $  SFPERATFO  FROM  1NCI0FNT  BAR  -  ANALVS'C 
1  TERMINATFD) 

986  FORMAT  ( 1H0.65HSPECIVFN  HAS  SFPFRA7F0  PROM  TRANSMITTER  BAR  -  ANALY 
ISIS  TERMinatFO) 

150  CONTINUE 

READ  (5.2)  XX< 

IF  (ICKK)  200.205.210 
999  STOP  „  . 

END  B"4 
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MRFTC  STORF  LIST 

SUBROUT I NF  STORE 

SUBROUT  I  NT  STORE  CONVERTS  THF  NON-D I  •-•F^S  I  O'.AL  QUANTITIES  INTO  THE 
SAME  UNITS  U5FP  FOR  INPUT.  THIS  CONVERSION  MAY  BE  SUPPRESSED  BY 
ENTERING  A  1  IN  COL  15.  CARD  7.  THF  OUTPUT  WILL  THEN  BE  IN 
THE  DIMENSIONLESS  FORM  INDICATED  IN  THE  PPOGRAM  DICTIONARY. 
REFLECTED  AND  TRANSMITTED  STRAIN  VALUES  ARE  ALWAYS  IN  ThE  ACTUAL 
UNITS  OF  STRAIN  { IN/IN). 

COMM-  .  A,A1 »A2.A3,AA,A5,ARrAl  .ARC  A?, o, PI  ,0  2, P3tn4,c5,RF, CI2D1 ,1 2) . 
1CI .C2.CBAR 1 .CBAR2 .CF.CO (1 1 »  50  ) .'"R.CS.rw.Ol »02»DC»DS,DT»DX,E (201  ,J2 
21 .El  «E?,FF(12) .F0< 11 ,50) ,EPS1 I2oo» ,EPs?( 1 2 1 .FP.EROt  50) .ES.ET  1 1?  )  , 
3ET0I50) ,EW,G. I . 1 1 . I  0 ( 12 ) .J.K.L .LL.LLL  < 1 1 i , L*ax , LMAX 1 , LMA X2 , « . N , NN , 
4NMAX ,NPAG£  »N STORE  »P I ?oi ,12) ,P1  ,P2.P0n 1 ,5ft) , PP » PR , PW .RH01 ,RH02, 
5Rh0S . S ( 201 .12) . SLOPE, SMAXI 201  ) ,S0(11  .SO) .SR.SS.SW.sy.T ( 20) ,12) »T1» 
6T2  ,  TIM (200) ,T0( 11 ,50) .TR.TS.UI  201  ,22) ,U0( 1 1 ,50) .UP > US, XBAR( 11 ) , XO 
<N  *  N-2 
DO  605  N*1 ,xr 
NSTORF  «  NSTORE  ♦  1 
DO  600  LOUT  *  1 , I 
L  «  LLL(LOUT) 

IF  ( M  .EO.  15  GO  TO  902 
T0( LOUT ,N$70RE )  =  T IL  ,N  *  *  A3 
SOI  LOUT .NSTORE )  =  SIL.N5  *  SY 
POILOUT, NSTORF)  =  P(l.N)  *  A2 
UOI LOUT, NSTORF)  =  UIL.N)  «  A2 
EOILOUT .NSTORE )  =  E(L.N)  «  A1 
COILOUT, NSTORF )  =  C(L,NI  *  CS 
GO  TO  ROO 

902  TOILOUT, NSTORE)  =  T ( L  » N ) 

SOILOUT, NSTORE)  »  S(L.N) 

POILOUT, NSTORE !  *  P(L,N) 

UOlLOUT, NSTORE !  «  UIL.N) 

EO I  LOUT ,N  6T0RF )  =  ElL.N) 

COILOUT .NSTORE)  =  CIL.N) 

900  CONTI NUF 

FRO  I NSTORE  *  *  E^IN) 

ETO I NSTORE )  *  FT  I  N  ) 

905  CONTINUE 

IF  INSTORE  .EO.  50)  C-0  TO  901 

PETU9N 

ENTRY  OUTPUT 
9P 1  NPAF.F  =  NPAf.F  +  i 
DO  920  LOUT  =  1,1 

WRITE  (6,910)  XBARILOUT),  NP AGE 

WRITE  (6.911)  I TOILOUT, JJ j , SO  I  LOUT , JJ ) .POILOUT, JJ) ,UOl LOUT, JJ) , 

1  EOILOUT.JJ) .COILOUT, JJ) ,  JJ=1, NSTORE) 

920  CONTINUE 

WRITE  (6.930)  NPAGE 

WRITE  (6.931)  I  TO! 1 ,JJ) ,ERO( JJ) ,  JJ= 1 , NSTORF) 

LOUT  •  L L L  1 1 ) 

WRITE  (6.940)  NPAGE 

WRITE  (6.931)  ( T0( I , JJ) ,ETO( JJ ) ,  JJ=1,NST0RD 

NSTORE  «  0 

RETJRN 

910  FORMAT  (  1  H 1  ,  6HXBAR  =  ,  F  7 . 5  ,«9X  ,4HP  AGE  ,  I  1  HO  ,  1 1  X  ,  4H7  I  ME  .  1  6x  ,  6HSTRES 
1S.1 5X.3HPHI , 1  OX, 1 7HPART ICLE  VELOC I T Y , 9X , 6HSTRA I N , 1 2 X , 1 OHWAVF  SPEFD 
2) 

911  FORMAT  ( 6E20.8  ) 

930  FORMAT  (  1  HI  ,  1  2HI  NF  I  DFNT  0  49 , 1  OPX  ,  4HP  Ar.r  ,  j  4,  ,  j  wp  ,  )  1  x  ,  4HT  I  Mr  ,  1  OX  , 

\  1 6HRFFL FCTFD  STRAIN) 

931  FORMAT  (2E20.8) 

940  FORMAT  ( 1H1 .15HTRANSMITTER  BAR ,9  7X , 4rtPAGF .14/1H0.11X.4HTI ME , 9X  » 

1  18HTRANSMITTEP  STRAIN) 

END 
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* IRFTC  ROUND  L  f.T 

SUBROUTINE  BOUND 

SUBROUTINE  BOUND  SOLVES  THE  CHARAC TTR : ST  IC  EOMATIONS  AT  THE 
BOUNDAR  I F  S  OF  THF  SPECIMEN.  THIS  SUOROL'T  I NF  ASSUMES  THAT  THE 
SPECIMEN-PRESSURE  BA J  INTERFACE  Is  ERJCTIONLFSS  AND  THAT  THE 
AXIAL  FORCE  AND  VELOCITY  OF  THE  SPECIMEN  AND  PRESSURE  BAR  ARE 
FOUAL. 

COMMON  A,A1  ,A2,A3,A4,AS,ARFAl  ,arca2,r,o1  ,f>2,R3,R4«PS,BF,r(?r>l  ,  12)  . 
1C1 »C2.CBAR1 ,CBAR2  «CE  *CO( 1 i . 5 0 ) ,CR ,C S ,CW ,0} ,02 .DC , DS ,DT ,DX , E ( 201 , 1  2 
2) »E1 ,E2»EF( 12) ,C0(11 .50) .FPS1 I ?oo ) »EPC2( 12) .ER.EROI 50) .ES.ET  1 12 )  . 
3ETOI50) ,EW,G» I . 11 .)D{ 12  > ,J.K.L,LL.LLL( 11 ) , LMAX , L MAX  1 , LMA X2 . M .N , NN » 
ANMAX.NPAGE .NSTORE .P( 20!  .12 )  .PI .P?.o0( I  1  .50) , PP , PR , PW.RH01 .RH02. 
5RH0S . S ( 20 1 .12) . SLOPE. SMAX  <  20 1  I . SO (I  1 ,50 ) , SR . SS . SW , SY . T ( 201 , 1  2 ) . T 1  , 
6T2. Tl M (200}, T0(!1 .50), TR.TS. 0(201,12) .00 1  11,50) ,UR . US, XBAR (1 1  I ,X0 
IF  (L-l)  700,700, 8op 
7 dp  T!1»N)=T( 2 ,N— 1 )  *  DX/CI2.N-1) 

IF  (T< 1 ,N)-C2*TIM( Jl)  7PA.70A.7O2 
70?  J=J+1 

SL0PE  =  C1*(EPS1 ( J)-FPS1 ( J-l ) )/(C2k  t  T!M( J)-T IMf J-l  )) ) 

70A  EPS2IN)  =  C1*EPS1  i  J-1I  +  SL0PE*  (  T  (  1  ,M-C2*T!M(  j-i  )  ) 

705  NN=N-1 

P1=RH0S*AREA1/(Rh01«CPAR1 ) 

P2  =P ( 2  »N- 1 ) — U ( 2 , N-l  )  +  2.*CBADl*FPS2(N)*FS/SY 

CALL  rrop 

P( 1 ,N I =PP 

PW«PP 

CALL  WAVE 

S(1.N5«SW 

CI1 cN ) =CW 

E(1 ,N)=FW 

U< ! ,Nja2.»CBARl*FPS2(Nl«FS/SY-P1*SW 

FF (N)*0.5*(S( 1 .N)*ARFtl *SY/E]-U( 1 ,N ) * SY/ [ CBAR 1 #F$ , j 
GO  TO  850 

BOO  TILMAX.N)=T (LHAX-; , N-l )+DX/C(LMAX-l ,N-1 ) 
d1=RHOS»ARPA2/  (°HC>2*FnAR2  ! 

D2*U(LMAX-1 «N— 1 ) *P( LMAX-1 ,N-1 > 

CALL  PROP 
PW«PP 
CALL  WAVE 
P ( LVAX  »N J »PW 
S I L  MAX , N ) 

C(LMAX.N)»CW 
E(LMAX,N)=EW 
U(LMAX,N)»P1*SW 
ET ( N) «S( LMA X .N)*A PE A 2* SY/F2 
850  CONTINUE 
RETURN 
END 
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MBFTC  STRESS  LIST 

SUBROU1  I  NT  STRESS 

SUBROUTINE  STRESS  SOL  VI  s  jmc  C  h  t.  PA<"  TrR  •  c  t  >  r  EQUATIONS  AT  ALL 
INTERIOR  POINTS  OF  T“F  cprop/pN  A^stOM  Nr  NO  INITIAL  OmANGE  IN 
CROSS-SrCTJ ON A i  ART a. 

OO^VON  A,Al,A2,A3,AA,A5,ARrAl,ARrA2,n,r-l,Q2,P;,OA,nS,OF,f(?f>l,l?), 

1C1 .C2.CBAR1  .CRAR2  »Cf  .COt  1 1  . 5o  )  ,CP  ,F S ,C  v  ,01  ,Q2  .F>C  . OS  ,r>T  ,f)x  ,e  ( 2o  1  ,  1 2 
2)  .El  .F2.FFI  12  )  .fOdi  .SO)  .rPa]  (  2001  .EPS2I  12  )  ♦  FR  *  FRO  (  50)  .Fo.ET  (  1  2  1  ♦ 
3ET01  50) .EW.G. I . 11 . ID! 12) .J.K.L.LL.LLL I  1  1  )  ,  LMAX  ,L«AX  1  , L  ,/AX2  »  v  *N  »  NN  » 
4NMAX.NPACE.NST0RF .PI201  ,  1 2  >  .  P 1  , P2 , PO Ill . G0 ) ,PP ,PR ,Pv.  ,PH01  ,RH02, 
5RH0S ♦ S ( 201 .12 ) . SLOPE tSMAX <  20 1  ) . SO  I  1 1 . 50 ) .SP . SS . S» . ST . T I  201 . 1 2 > . T 1 , 
6T2.TIMI20O)  ,T0(  1  1  .50)  ,  TP,  IS.  U(  201 ,12)  ,-jO<  1  1  .  50  )  .UR  ,US,  Xf»  AR  (  11  )  ,  XO 
T1=T(L-1.NN  M-IOx/OL-l ,NN  )) 

T2*T ( L+l «NN  )+(OX/C(t +1 ,NN  1) 

IF  ( A8S  (  T 1  /  T2  -  1.)  .LF.  o.aoooaoI)  GO  TO  420 
IF  ( T 1  —  T  2 )  A0O.A20.A1O 
AOO  T ( L  »N )  =  T  1 
LL=L+1 
CALL  SHIFT 
GO  TO  A 1 S 
A10  T ( L  »N ) =T2 
LL  *L-1 
CALL  shi^t 

IF  ( L-? )  A15.A11.A15 
All  IF  ( TR-C2*T I M ( J-l ) )  A12.A13.A13 
A12  J  =  J-1 

SLOPE  =  C1»(FPS1 ( J)-EPS1 ( J-l ) )/(C2«I  TIM,  D-TIMIJ-1  )  )  ) 

A 1  3  FPS?(N)=ri*FPSl  I J-l )+n  *(FPS1 ( Jl-FPSl ( j-l ) )•  I TR-C2*TI'i|  J-l  )  )/ 

1  (  C2*  (  T  JM(  J  )-T  !**{  J-l  )  )  ) 

415  PILL  »NN ) *PP 
SILL ,NN)*SP 
E I LL ,NN ) «ER 
UILL.NN)=UR 
C I LL  »NN ) sCR 
T  5  LL  *NN ) *TP 
GO  TO  A50 
420  TIL.NJsTI 
A50  CONTI NUF 

IF  IL.E0.1I  GO  TO  Aeo 

PIL.N)=.5*IPIL+1 »NN  J+PIL-I.NN  )-U(L+l.NN  )+UIL-l.NN  )) 
U(L»N)«.5*(UtL+l»NN  )+UIl-l,NN  1-PtL+l.NN  )*?IL-1»NN  )) 

IF  IABSIPIL,N)/P(L,N-1)  -  1.)  .LF.  0*0000001)  GO  TO  47R 
IF  IPIL.NJ-PIL ,N-1 ) )  A55.A7s.A75 
A 55  S(L,N)=0.5*|SIL+1 .NNM-s iL-1 ,NN)+UI L-l .NN)-U«L*1 *NN1 ) 

UIL  .N) rn.5*I SI L-l ,NN)-S(L  +  1 ,NN)+UIL-1 ,N\  ) ♦ U I L ♦ 1 ,N*M  ) 

A70  PIL»N)*SIL,N) 

CIL.N)*1. 

EIL»N)*FlL»N-l ) ♦S I L ,N )— S I L  »N— 1 ) 

GO  TO  A 60 
A 75  PW*P I L ,N ) 

CALL  WAVE 
sil»m*sw 

C I L  »N ) *CW 
E I L  »N ) =Ew 
return 

A 79  SIL.N)  =  SIL.N-1 ) 

EIL.N)  =  ElL.N-11 
CIL.Nj  *  CtL.N-l) 

480  CONTINUE 
RETURN 
FNO 
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APPENDIX  B 

SIPFTC  WAVri  LIST 

SUBROUTINE  WAVE 

SUPROUT 1 NE  WAVE  IS  THE  ONLY  SUBROUTINE  WHICH  REOUIPES  USE  OF  THE 
ST°F  SS-STP A  I N  RELATION  0r  THE  SPECl^N.  T“T  RFLAT  !  ONSHI  P  USED  IN 
THF  PLASTIC  PEG  I  ON  IS...  r  s  PI  P?»S  +  B3*S**2 
THIS  RELATIONSHIP  HAS  ALREADY  BEEN  NORMALIZED  SO  THA T . . B 1 +B2+B3= 1 
THIS  ROUTINE  CONSISTS  OF  THREE  SEPARATE  PARTS. 

COMMON  A.A1  ,A2.A3.A4,A5,APFA1 .ARRA2.0.R1  » n 2 . B 3 , RA  ,  p 5 , BF , c ( 2 0 1 » 1 2 ) . 
IC1  ,C2,CBAR1 .C8TR2.CE ,CO(ll .50) »CR .C S .f w . D1 .02 .DC . DS .OT ,Dx ,E t 201 . 12 
2) .E1»E2.FF( 12) ,F0( H .SOI .EPSI l  2op ) «EPs2t 12) .ER.EROl 50) .ES.ET 112  I . 
3ETO(50)»EW*C-,I,I1»ID(12)«J.K.L,LL.LLI  (11)  »  LMAX  .  LMA  X 1  . LMAX2 ,N , NN , 
4NMAX.NPAGE.NST0PE ,P( 2M  ,12)  ,P1  ,P2  .POn  1  .so  )  .PP.PR  .Pw.RHOl  »Ph02* 
SP^OS.Sf  701,1?)  .si.OPR.sma*  I  ?01  )  .SOI  11  .SO)  , SR ,SS, SW . SY ,T t 20]  .1  2)  ,T1  . 
6  T  2 . T I M ( 2  O  s ) ,io< 11  ,50) .TR,T$»U(201 ,12) «  UO (11.50) ,UP » US. XftAR ( 1 1 ) ,X0 
FOR  A  GIVEN  IMPACT  VELOCITY  THIS  PART  OF  THE  SUBROUTINE  WILL 
COMPUTE  THE  SPECIMEN  STRESS.  STRAIN,  AND  WAVE  SPEED. 

IF  < (PW/SMAXIL) )  .GT.  0.9999999)  GO  TO  210 
200  SW=PW 
CW=1. 

EW  =  r ( L .N— 1 )  +  SW-S 1 1  . N- 1 ) 

GO  TO  2P0 

210  SW=( ( I 3.*B3*(PW-1 . ) )+(P2+2.*P3)**l.S)«*(2./3.)-P2!/(2.*P3) 

CW=1 . / SORT (P2+2.*P3*SW) 

EW  =  B1-»B2*SW+B3»SW**2 
200  CONTINUE 
RETURN 

FOP  GIVEN  CHARACTERISTIC  VALUFS  AT  THE  BOUNDARIES.  THIS  PORTION 
OF  THF  SURROUT I Nr  WILL  COMPUTF  THE  STPFss  AND  IMPACT  VFLOCITV. 
ENTRY  PROP 
IE  ( L-l )  301, SOI, 302 

3nl  IF  {(P2/AA)  .C,T.  o.oooooao)  GO  TO  3o5 

303  P2  =P2— P  ( 2 »N- 1 >  +  S I  2  »N— 1  ) 

GO  TO  300 

302  IF  ( { P 2 / A 5 )  .GT.  0.9000001,  GO  TO  306 

304  P2=P2-p(LvAX-l , N— 1 >+SCLMAX-i.N-l > 

GO  TO  300 

SOS  AA=P2 

GO  TO  ’10 
306  AS  =  P2 

GO  TO  310 
300  SP  =  P2 / I  PI  +  1 .  ) 

PP»SP 
GO  TO  300 

“10  V=P 2-1 .♦?“?♦?. *«A}**i ,s/( 0.*OS) 

X=1.S*P2/°y—  C°.*P1**?)/(S.,P0» 

Y=.75*<B2/83)*«2+(9.*Pl*V)/(4.*B3t 

Z=.12S*(B?/B3J**,-(0.*v*V5/(8.*B3) 

SP*R0CT(X,Y,Z) 

3S0  PP=l.+t(B2+2.*B3*SP)**I.S-{B2+2.«B3)**1.5)/(3.*B3) 

300  CONTINUE 
RETURN 

FOP  A  GIVEN  Tlwr  pasTO  ON  IMF  SLOPE  OF  CNF  CHARACTERISTIC.  THIS 
PORTION  OF  THr  SUBROUTINE  „ILL  COMPETE  THE  VALUE  OF  THE  OTHER 
oEojiPED  characteristic, 
entry  SHIFT 

DT*T(LL.NNI-t«LL,NN-1  ) 

IF  (  (  SILL.NN-l  ) /SMAXI  LL  ))  .GT.  0.9999099)  GO  TO  MO 
900  JF  (  (S(LL,NN)/SVAX(LL)  )  «  r-  7  »  0.999O9OO)  GO  TO  402 
401  TR«TCL.N)-DX 

po«R( ll ,nn-i )♦{ tr-tu l ,nn-i ))*; pill, nn) -pill »nn-i ) > /dt 
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APPENDIX  B 

UR=U(LL,NN-1  )  ■*■  (  te-t  U.  L  »  NN-1  )  )  *  {  U  (  LL  .  NN  ) -U  (  LL  ,NN-J  )  ) /DT 
GO  TO  480 

402  TR  =  T  <  LL  •  NO- 1  )+DT*  (  SVAXC LL  >-P( LL.NN-1  J  )  /  (  P  (  LL  »  NN  )  -P  t  LL  .  NN-1  )  ) 

CSMAX  si./  SORT ( 82  +  2. *D3«SMAX< LL ) J 

PSMAX  =  1 .*( ( 82+2.*«3*SMAx(LL ) )»»1 ,5-(R2+2.*R3)«*l .5)/ ( 3.*B3) 

IF  (TR  +  DX-TU.N)  )  403.401,401 

403  IF  (TR+OX/C^^AX-T f  L  »N  ) )  406,404,404 

406  X  =  0. 

Ys(DX/t  TILL,  NN1-TP)  )*(1./C  (LL  ,  NN )  **  3- 1  . /CS'-'A  X*  «  3  ) 

Z*-l./CSKAX«*3+<  (TR-T(L.N)  )/{T(LL,NN)-T0>)*(1./(C(LL,N\)«*3I- 
1  l./CSMAX**3) 

CRI  =  ROOT { X , Y . 7 ) 

TR1  =  T(L.N)  -  DX  *  CPI 

PR  =  PSMAX  +  (  (  T  R  1  -  TRJMTU  t  ,\N1  -  TP  t  I  «  (  P  (  L  L  ,  N  N )  -  PSMAX) 

TR  =  TP1 

OR  =  'JILL,  NN-1  >♦{ TR-T (LL ,NN-1 >)« (U(LL .N')-U(LL .NN-1 )) /DT 
GO  TO  480 

404  PRsPSNAX 

ORsUC  LL  .NN-1  1  ♦!  TD-T  (LL  ,NN~i  ))  *  (  U(  LL  ,NN  )-u(  LL  .NN-;  )  )  /DT 
GO  TO  480 

410  IF  ( (SiLL.NN)/Sv4X(LL > )  .OT.  0.09990°°)  GO  TO  4?A 
415  TR  =  T(1 ,N)  -  OX/CILL ,NN-1 ) 

S°  —  c ( L  L , NN— 1  ) 

°R  =  P 1 1 L .NN-1  ) 

U°  =  U(LL.NN-I) 

CR  *  CILL.NN-1 ) 

ER  =  E I LL.MN-l) 

GO  TO  495 
U?-  X=0. 

Ys (DX/OT )*(1 ,/C (LL  »NN)*«3-1  ,/C (LL ,NN-1 )*»3) 

2  =  -l  ./CUL.NN-l >**3+( ( T (LL .NX-1 ) -T ( L ,N ) ) /DT ) * ( 1 . / ( C ( LL ,NN ) **3 )  ./ 

1 (CILL.NN-1 )**3  )  ) 

CRI sROOT ( X , Y , 7 ) 

TR  =  T(L,N!  -  DX  «  CRI 

PR  =  PUL, NN-1)  ♦  (iT»  -  T  ( LL  ,NN- 1  )) /DT  )  *  (P(LL.NN)  -  PlLL.NN-1)) 
UR  =1*  <  LL ,NN-1 )«-(  TP-T  (LL.NN-1 ) i * ( U C LL ,NN )-U ( LL .NN- 1 ) ) /DT 
4S0  IF  (  (PR/SF"AX(LL))  . GT .  0.°9999°9)  GO  TO  490 
485  SR=PR 
CR  =  1. 

FRst'UI.  .NN-1  )-*-S°-'(LI  , NN—  i  J 
GO  TO  495 

400  SRs((3.*fi3*(PR-l.)+(P2^?.*P3)**l.c)«* ( 2./3. )-°2 > / I 2.*B3 ) 

CP*'  1. /SORT  (B2*2. •F>3*°=) 

ER=B1+P2*SR+D3*SR**2 
495  CONTINUE 

IF  (LL.fO.l )  GO  TO  498 
IF  (LL.rO.LMAX)  AO  TO  499 
DETURN 

408  FE  (  NN  )  =EF  (  NN-1  l  +  (  T°-T  (LL  ,NN-1  )  )  •  (  FF  (  NN  )-FF  (  NN-1  )  l  /r>T 
RETURN 

4  9°  ET  (NN)  =  ET  (NN-1  )  +  (  (R-T  (  LL  ,*.N-1  )  )  *  (  E  T  (  NN  )  -t  T  (NN-1  )  )  /DT 

return 

END 
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$  I PFTC  POOT  LIST 

FUNCTION  ROOT  IX,Y,2> 

FUNCTION  ROOT  FINOS  THE  ONE  REAL  ROOT  OF  THE  CUBIC  EQUATION 
X»*3  *•  A*X**2  +  8* X  +  C  =  0.  THE  ARGUMENT  LIST  OF  THE  CALL 
STATEMENT  CONTAINS  THF  THREE  COrFFIClFNTS  A,  P,  AND  C. 
DOUBLE  PRECISION  A,B,0,O,F,F 
A=Y/^.-( X/3. )  *»2 
8  =  <  X  / 3.  )»*A  — x*Y/6.4-7/2. 

C=B**2+A** a 
IF  (C)  10*20.20 
10  PH I=-8/0SQRTI-A**3) 

ROOT«-X/3.+2.*DSORT(-A)*COS{PH!/3. ) 

GO  TO  30 
2o  0*— P+OSORT { C ) 

F«3+DS0RT(G  ) 

F=1 ./3. 

IF  (D)  25,22*22 

22  IF  IFI  2  3  »  24  »  2 A 

23  RGOT*-X  /  3.+D**F-*{-E)**F 
GO  TO  30 

24  POOT*-X /3 • +0**F-E»*F 
GO  TO  30 

25  ROOT =-X/3.-(-D)**F-E**F 
30  RETURN 

END 
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APPENDIX  C 

CALIBRATION  OP  RECORDING  DEVICES 

The  measurement  of  stress  and  particle  velocity  in  each 
pressure  bar  is  accomplished  indirectly  by  recording  transient 
electrical  signals  generated  by  the  change  in  resistance  of 
etched  foil  strain  gages.  This  technique  differs  somewhat  from 
the  methods  employed  in  the  static  measurement  of  strain  where 
the  voltage  required  to  "balance"  a  Wheatstone  bridge  is  propor¬ 
tional  to  the  strain.  The  transient  measurement  of  strain 
requires  the  bridge  to  be  initially  balanced  and  the  amount  of 
"imbalance"  recorded  as  the  strain  is  applied. 

Each  strain-measuring  station  consists  of  two  strain  gages 
which  are  connected  to  opposite  arms  of  a  Wheatstone  bridge  as 
illustrated  schematically  in  figure  C-l(A).  Each  pair  of  strain 
gages  is  connected  to  a  separate  channel  of  a  transducer  input 
conditioner  which  contains  six  Wheatstone  bridges  (B  and  F  Model 
1-220B4).  The  system  is  calibrated  prior  to  each  test  by  switching 
known  precision  resistances  into  the  bridge  circuits  and  recording 
the  resulting  deflections  on  Tektronix  Model  531  oscilloscopes 
equipped  with  TyPe  D  plug-in  units.  The  intensity  is  modulated 
every  50  microseconds  by  means  of  a  separate  external  oscillator 
which  permits  calibration  of  both  the  vertical  and  horizontal 
axes  of  the  oscilloscope  as  shown  in  figure  C-l(B). 

Hie  equation  which  describes  the  output  of  a  bridge  with  four 
initially  equal  arms  and  two  active  gages  in  opposite  arms  is 
given  by: 


AV  i  AR 

*r  =  *  ?r 


(1 


n) 


(C-l 


> 


where  the  (1  -  n)  factor  is  required  to  account  for  tne  non- 
linearities  which  arise  for  large  changes  in  resistance.  This 
correction  factor  is  given  by: 

(1  -  n)  =  (1  +  (C-2) 

The  "engineering"  strain  is  proportional  to  the  ratio  of  the 
change  in  resistance  to  the  original  resistance  of  the  gage.  The 
constant  of  proportionality  is  the  gage  factor,  G,  supplied  by 
the  manufacturer.  Equations  (C-l)  and  (C-2)  may  be  rewritten  in 
terms  of  strain  as: 

%  =  *Gc  (1  -  n)  (C-3) 
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ana 


(1  -  n)  =  (1  +  Ge)'1 


(C-*) 


Since  the  gage  factor  has  a  nominal  value  of  two,  straina  on  the 
order  of  2500  microstrain  may  be  measured  to  within  £  percent  when 
the  (1  -  n)  factor  is  neglected.  Since  the  strains  in  the  pressure 
bars  never  exceeded  this  value,  n  was  assumed  to  be  equal  to  zero. 


The  calibration  resistors  are  wired  in  parallel  with  only  one 
of  the  active  arms  of  the  bridge,  figure  C-1(a).  Therefore,  the 
change  in  resistance  of  the  bridge  is  only  half  that  which  will 
occur  when  both  gages  are  strained.  The  equivalent  strain 
corresponding  to  a  given  calibration  resistor  is  given  by: 


cal 


* 


cal' 


Prom  equation  ( C— 3 )  we  obtain: 


(C-5) 
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^cal 


\  0 


cax 


(C-6) 


Therefore,  the  strain  is  related  to  the  measured  change  in  voltage 
by  eliminating  V  and  G  from  equations  (C-3)  and  (C-6)  to  obtain: 
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Each  strain  corresponding  to  a  calibration  resistor  may  be 
computed  from  equation  (C-5)  and  the  corresponding  displacement 
on  an  oscilloscope  determined  experimentally  as  shown  in  figure 
C-l(B). 

The  elastic  wave  speed  in  each  pressure  bar  is  determined  by 
Impacting  each  bar  separately  and  observing  the  multiple  reflections 
of  the  propagating  wave  on  an  oscilloscope  sweeping  at  one  milli¬ 
second/centimeter.  The  output  of  a  calibrated  external  oscillator 
is  superimposed  to  obtain  an  accurate  time  base.  An  accurate 
measurement  of  the  length,  A,  of  each  pressure  bar  and  the  total 
elapsed  time.  At,  between  (n  +  1)  peaks  observed  in  the  propagating 
wave  permits  the  average  wave  speed  to  be  computed  from  the 
expression: 


„  .  2  n  tQ 

W  —  11.  1  .  ■ 

0  At 


(C~8) 
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speed  for  the  titanium  pressure 
was  found  by  this  method  to  be 
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(A)  SCHEMATIC  DIAGRAM  OF  TYPICAL  CHANNEL  OF  TRANSDUCER 
INPUT  CONDITIONER 
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(B)  TYPICAL  CALIBRATION  RECORD  (P1  =  ©6.7  x  10'6  IN/IN; 
?2  =  817.5  x  10-6  IN/IN;  P3  =  2366.7  x  ICf6  IN/IN) 


FIG.C-1  CALIBRATION  OF  STRAIN  MEASURING  EQUIPMENT 
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